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Abstract 


Currently there is no definitive description for the accelerated ex¬ 
pansion of the Universe at both early and late times; we know these 
two periods as the epochs of inflation and dark energy. Contained 
within this Thesis are two studies of inflation and one in the con¬ 
text of dark energy. The first study involves two noncanonical ki¬ 
netic terms each in a two-field scenario, and their effects on the 
generation of isocurvature modes. As a result, these terms affect 
the isocurvature perturbations produced, and consequently the Cos¬ 
mic Microwave Background. In the following study, the impact of a 
sharp transition upon the effective Planck mass is considered in both 
a single-field and two-field model. A feature in the primordial power 
spectrum arising from these transitions is found in single-field mod¬ 
els, but not for two-field models. The final model discussed is on the 
subject of dark energy. A type of nonconformal coupling is exam¬ 
ined namely the “disformal” coupling; in this scenario a scalar field 
is disformally coupled to matter species. Two consistency checks are 
undertaken, the first to provide a fluid description and the second, 
a kinetic theory. From this, observables are constructed and used 
to create constraints on the individual coupling strengths. 
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Chapter 1 


Introduction 


Cosmological inflation has been widely accepted as the theory giving rise to our 
Universe and over the decades, cosmologists have attempted to understand the 
mechanism behind this. More recently, evidence of a second period of acceler¬ 
ated expansion in our Universe has been found. However, the reasons behind 
this late-time acceleration are unknown. It could be as a result of a new energy 
form, namely dark energy, or clue to modifications of gravity. 

This introductory chapter is presented as follows: first, a description of the 
three tests of General Relativity which is proceeded by a summary of the dis¬ 
coveries made in the field of cosmology. Following this, the problems in cosmol¬ 
ogy will be presented alongside its solution — inflation. Lastly, a description of 
cosmological perturbation theory is given before ending with an outline of the 
Thesis. 
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CHAPTER 1. INTRODUCTION 


1.1 General Relativity and the three tests 

The Theory of General Relativity by Einstein was a culmination of work at¬ 
tained during the period from 1905 to 1915. Many papers detailing the technical 
nature were published during this time, but it was not fully presented until 1915 
and 1916 [HI [9]. However, tests were required to scrutinise his theory and three 
were proposed: the perihelion precession of the planet Mercury, the deflection of 
light by massive bodies and gravitational redshift of light near compact massive 
bodies. 

We will work with the mathematical framework provided by the Theory of Gen¬ 
eral Relativity, the most successful theory describing gravity in curved space- 
time. In order to study models which are extensions of General Relativity, it 
must be ensured that these models satisfy the three historical tests. 

The first test — Mercury’s perihelion precession 

During the 19th Century, a French mathematician named Le Verrier undertook 
an extensive study into the nature of the planets. Through the usage of New¬ 
tonian mechanics, Le Verrier had discovered that when all the motion of that 
surrounding planets had been accounted for, it was not enough to explain the 
perihelion precession of the planet Mercury. The unaccounted precession was 
calculated to be approximately 38 arcseconds per century [10J. In his Letter 
of 1859, he gave possible reasons for this discrepancy which included an undis¬ 
covered planet and objects, such as asteroids, orbiting between the Sun and 
Mercury. However it was not until November of 1915 that Mercury’s missing 
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perihelion precession would be solved by the Theory of Relativity. In 1915, Ein¬ 
stein had calculated it to be 43 arcseconds per century, which was in agreement 
with the current-day observational values of 45" ± 5" per century m- 

The second test — Light deflection by massive bodies 

Another test involved the deflection of light in the presence of a massive object. 
This massive object induces the light deflection by causing a change in the 
curvature of space-time in its locality. The bending of light near a celestial 
object, in this case the Sun, was first stated and shown by Soldner in 1804 
through a calculation using Newtonian mechanics m Through his calculation, 
Soldner found the amount of deflection of the light ray caused by the Sun is 
0.84”, and this figure was also calculated by Einstein in 1911 H3J by using 
the equivalence principle. The equivalence principle is the assumption that the 
gravitational force experienced by a local observer is the same as that of an 
accelerating reference frame. However, Einstein corrected the calculation in 
1916 by using his Theory of Relativity and from this, found that the amount 
of light deflection doubled to 1.75”. This led to the solar eclipse experiments of 
May 29, 1919 performed by Eddington and Crommelin [T4J whose purpose was 
to confirm if gravitation had an effect on light and if yes, whether Newtonian 
mechanics or Einstein’s Theory of Relativity was the correct description of 
gravity. It was noted that this eclipse would be particularly favourable due to 
the unexpected number of bright stars in the background, which gave many test 
stars to compare against. The path of the eclipse was to pass over Northern 
Brazil and the Atlantic Ocean, eventually reaching Sao Tome and Principe, an 
island country off the West Coast of Africa, before passing over the continent. 
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As a result, two expeditions were commissioned: one to be based in Sobral in 
Northern Brazil led by Crommelin, and the other on the island of Principe to be 
led by Eddington. Both experiments performed were found to be in agreement 
with the prediction made by Einstein’s Theory of Relativity, with results from 
the experiment led by Crommelin found to be 1.98" ± 0.16" and Eddington 
1.61" ± 0.4". 

The third test — Gravitational redshift of light 

The gravitational redshift of light close to compact massive bodies was the last 
test proposed by Einstein. It was confirmed by Adams in 1925 mm through 
the study of a binary star system in Sirius. Previous to this, it was also Adams 
in 1915 who discovered that Sirius A had a faint companion, specifically, a white 
dwarf - the first to be found, and was named Sirius B im. It was Eddington 
in 1924 [18] that proposed that these white dwarfs are of very high density in 
comparison to other stellar objects; his statement was made whilst specifically 
discussing the binary system of Sirius. Both the effective temperature and ab¬ 
solute magnitude of Sirius was used to calculate the radius, and in combination 
with an assumption of a mass range, came to the conclusion that white dwarfs 
are extremely high in density. Alongside this discussion, Eddington also gave 
an alternative explanation in which the stellar object has instead a very low 
effective temperature, and the ability to generate the same spectrum which an 
observer must be satisfied with. In order to determine one of these two ideas, 
the gravitational redshift of the spectrum generated by the white dwarf must 
be measured. For Sirius B, a very high density stellar object, Eddington cal¬ 
culated the corresponding Doppler shift to be 20kms~ 1 . In the 1925 paper by 
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Adams, he had used a 100 inch reflector telescope - an upgrade from the 60 inch 
telescope used to detect Sirius B, to make spectroscopic observations. These 
observations and taking into account the radial velocity of Sirius led to Adams 
calculating the displacement of the spectrum of Sirius B to be 21kms -1 — in 
agreement with Eddington. 

There are other tests of General Relativity aside from the three historical tests 
stated; one such example uses binary pulsars, with the first system discovered 
in 1974 by Hulse and Taylor [19]. Observations of binary star systems have 
shown that the pulsar’s orbit contracts over time as a result of the emission of 
energy through gravitational waves. 

In 1917, Einstein’s famous paper titled “Kosmologische Betrachtungen zur allge- 
meinen Relativitatstheorie” [20] was published where he states the field equa¬ 
tions and from this, obtains a model where the Universe is both closed and 
static, through the term which he introduced as the cosmological constant. A 
positive cosmological constant was inserted into the field equations as a counter 
term to the effects of gravity — in a Newtonian model of the Universe, all the 
structures within will collapse under gravity, and the term being positive en¬ 
sured that the Universe was both closed and finite. Einstein and other physicists 
such as de Sitter [21] did not agree with the presence of the term. However, the 
term has been confirmed to be nonzero by supernovae observations performed 
in the late 1990s. 


We shall begin with a brief summary of the important discoveries made in 
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cosmology. 

1.2 Discoveries in Cosmology 

Expansion of space 

In Spring of 1917, an American astronomer named Slipher published a paper 
containing 4 years of study into the light emitted from surrounding galaxies ra- 
In this study starting in 1912, Slipher had stated that the size of the telescope 
is irrelevant in obtaining the spectroscopic measurements of the galaxies but 
instead, it was the camera of the spectrograph that had the greatest influence, 
specifically the lens speed. The spectrograph was attached to a 24 inch refractor 
telescope at the Lowell Observatory in Arizona, USA, and used to take between 
40 and 50 spectrograms of 25 spiral galaxies. By calculating the Doppler shift of 
the spectra lines, Slipher inferred the radial velocities of the galaxies and found 
that the average radial velocity was 570 kms -1 . In addition to this, most of the 
spiral galaxies were found to be moving away from the Solar System; Slipher 
had discovered that the light from the majority of the galaxies were shifted 
towards the infra-red section of the electromagnetic spectrum, known as red- 
shift. Hubble extended this study during the 1920s, resulting in his 1926 paper 
classifying the various types of galaxies in the Llniverse based on their shape, 
mass and luminosity [21|J. Further study into the various galaxies’ velocities and 
their distances led to the formation of the velocity-distance relation, presented 
in 1929 [23]. However, Hubble did not explicitly state the now-known relation 
v = H 0 d where v is the receding velocity, d is distance and H 0 is Hubble’s 
constant, and it was in 1927, a Belgian Roman Catholic priest named Lemaitre 
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who stated the relation between the recession velocity and distance [25] . and 
added the statement that it is natural for the galaxies to have a velocity which 
is receding with an expanding Universe. Furthermore in 1928, Robertson, an 
American physicist, also found a similar relation to Lemaitre [26]; Robertson 
presented a review of the important findings in relativistic cosmology in 1933 

( 23 • 


Two extensive studies into both the static and expanding Universe were first 
presented in 1922 and 1924, respectively by a Russian mathematician named 
Friedmann [28;, 29]. With these two papers, he proposed a space-time which 
is allowed to grow and/or shrink depending on time, and proceeded to write 
down the equations to which we now know as the Friedmann equations. The 
difference between the two papers is the sign associated with the curvature of 
the Universe; the first of the two papers, Friedmann considered a Universe with 
a positive constant curvature and hence, a static model with the Universe end¬ 
ing in a Big Crunch. In the other paper, a negative constant curvature leading 
to an expanding model. 

In the very same paper of 1927, Lemaitre independently obtained the same 
results as Friedmann. However, the work was written in French and published 
in an obscure journal in Brussels and hence, was relatively unknown to most 
of the scientific community; it was later published in English in 1931 [30] . The 
independent work performed by Lemaitre gave a greater appreciation for the 
two papers published by Friedmann. Lemaitre was searching for a compromise 
between the solutions given by de Sitter [31] and Einstein; the de Sitter so- 
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lution involved the assumption of no matter existing in the Universe and so 
the density equalling zero, whereas Einstein’s solution contained matter and 
led to a relation between the matter density and the radius of the Universe. 
In order to make a compromise between the two solutions, matter must exist 
in the Universe and adding to this, Lemaitre proposed that the radius of the 
Universe was allowed to vary. The solution he obtained in the paper revealed 
that the Universe naturally expands without bounds. Towards the end of his 
work, he stated that there still a need to find the reason for this expansion, 
which Lemaitre would do in 1931. 

The hot Big Bang 

The idea of the Universe starting initially from a single point — the Big Bang 
theory, was proposed by Lemaitre in his 1931 Letter [32] and the idea further 
developed in 1933 [33| [37 • 

One of the consequences of a hot and dense Big Bang is a remnant of a cooled 
thermal background radiation. The Cosmic Microwave Background (CMB) was 
discovered by Penzias and Wilson at Bell Laboratory in 1965, resulting in the 
award of the 1978 Nobel Prize for Physics [35]. A 20 foot horn reflector an¬ 
tenna was used to make the measurements and from this, a temperature of the 
background radiation was recorded to be 3.5 ± 1 K. The result was indepen¬ 
dent of direction revealing the Universe to be isotropic. The idea of a cooled 
thermal background was first postulated by Dicke and collaborators during the 
1960s which also led to a 1965 paper on the CMB temperature [36]; Penzias and 
Wilson cited this paper as it provided a possible explanation for the existence 
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of background radiation. Dicke and collaborators measured the temperature 
using a radiometer with a receiving horn attached which has the ability to 
measure radiation at the wavelength of 3cm. They found that with a wave¬ 
length of 3cm, the measured temperature of the cosmic blackbody radiation 
is 3.5 ± 0.5 K. In addition, the Letter was written with full acknowledgement 
of the work performed by Penzias and Wilson of Bell Telephone Laboratories. 
Both their works, led to the acceptance of the Big Bang theory. 


1980s — Inflation and new species (dark matter) 

There were three major problems of the Hot Big Bang model that surfaced by 
1980; they were namely the flatness, horizon and magnetic monopole problems. 
In order to solve the problems, the Universe must undergo a period of accel¬ 
erated expansion. Guth in 1981 outlined the problems clearly and proposed 
the idea of cosmological inflation [37J as a remedy for this; other authors who 
also independently worked on this concept are Albrecht and Steinhardt [38] and 
Linde [39] HO] , We will come back later to the theory of inflation. 

In the early 1930s Zwicky discovered that ordinary matter accounts for a very 
small proportion of the total matter in the Universe. By measuring the velocity 
dispersion of galaxy clusters through the viral theorem m, lie found that the 
visible matter was not enough to account for the fast movement of these galax¬ 
ies. He postulated that there must be some form of invisible matter contained 
within these clusters and this was later realised to be cold dark matter (CDM). 
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Several cosmologists in the 1980s including Peebles [42], Blumenthal [43] and 
Davis m, presented their findings which proposed that CDM was required for 
structure formation. Up to now, the models that were considered only con¬ 
tained ordinary matter (baryons, radiation and massless neutrinos) and were 
unable to produce the observed primordial fluctuation amplitude. The addition 
of the word “cold” was to affirm that the dark matter was nonrelativistic during 
the lifetime of the Universe. In the paper by Peebles in 1982 [42], CDM was 
included alongside ordinary matter and as a result, it caused the amplitude of 
the primordial fluctuations to change such that it could satisfy the current-day 
observations. 

Dawn of modern cosmological observations 

The Cosmic Background Explorer (COBE) in 1992 revealed that the back¬ 
ground radiation of the Universe is uniform on large angular scales. By the 
removal of various components such as the uniform background and dipole gen¬ 
erated by the Earth’s motion, primary temperature anisotropies were uncovered 
with an amplitude of AT/T ~ 10 -5 [45, 06]. The FIRAS instrument on COBE 
which investigated the blackbody nature of the CMB, produced the first precise 
measurement of the CMB temperature [47] recorded to be 2.726 ± 0.010 K; the 
most recent measurement was taken in 2009 [48] with 2.72548 ± 0.00057 K. 

In 1995, Goobar and Perlmutter, two members of the would-be Supernova 
Cosmology Project released a paper which discussed extending the apparent 
magnitude-redshift relation from z ~ 0.5 to z ~ 1 through the usage of Type 
la supernovae, and using this to measure the cosmological constant and the 
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density parameter f2 0 . Other papers followed which described the techniques 
and searches for the Type la supernovae between the redshifts z = 0.35 — 0.65 
j IT]. [50]. In 1998, observations of Type la supernovae showed that the cosmo¬ 
logical constant is nonzero and hence, the Universe is currently experiencing a 
period of accelerated expansion. This was found by two groups, the Supernova 
Cosmology Project [5TU52] , and Supernova Search Team [53, 54], independently 
in 1998. 

In 2001, the Wilkinson Microwave Anisotropy Probe (WMAP) satellite was 
launched for a mission of over 12 years, to take precise measurements of the 
CMB [55]; the final Nine-year WMAP data is presented in [55], Recently, 
Planck, a satellite launched by the European Space Agency, released their re¬ 
sults it found that the Universe is composed of 31.7% matter (of the total 
Universe, 4.9% is ordinary matter (baryonic) and 26.8% is cold dark matter) 
and 68.3% dark energy [57] , 

Alongside this observational evidence, the theories developed in the early 1980s 
by Peebles, Blumenthal et al. and Davis, this led to the establishment of the 
ACDM model; the concordance model in cosmology. A is known as the cosmo¬ 
logical constant representing dark energy in the Universe. It has been accepted 
that dark energy is the cause for the accelerated late-time expansion of the 
Universe. 

Up to now, there is no understanding of the mechanism that initiated inflation 
nor knowing the true nature of dark energy. Due to this, cosmologists have 


12 


CHAPTER 1. INTRODUCTION 


turned to alternative theories in order to explain both these predicaments. 


1.3 Preliminaries 

We shall use the mathematical formalism presented in General Relativity to 
describe the Universe. We know that the Universe is both homogenous and 
isotropic on large scales due to the observations of the CMB, and is expanding 
from Type la supernovae observations. This can be described in a mathemati¬ 
cally succinct form known as the Friedmann-Robertson-Walker (FRW) line ele¬ 
ment [58] I59J. The line element is maximally spatially symmetric and presented 
in the spherical coordinates. It is given by 



( 1 . 1 ) 


where a(t ) is the scale factor which is a function of the cosmic time t describ¬ 
ing the expansion of space, and k is the spatial curvature. The type of spatial 
curvature is assigned to a value of k; the values that k can take are k = — 1, 0,1 
describing an open, flat and closed spacetime, respectively. 

In general, a line element can be expressed by the following 


ds 2 = g^ v dx fl dx l/ , 


( 1 . 2 ) 


where g is the metric of the spacetime. In this Thesis, a flat space is used k = 
0, and with an expanding spacetime where the metric is g, iu = diag(—1, a 2 , a 2 , a 2 ), 
given in Cartesian coordinates. 
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In Einstein’s 1917 paper [20], a cosmological constant A was introduced in order 
to counter gravity and obtain a static Universe. The Einstein-Hilbert action for 
the case with a cosmological constant and additional matter fields Sm is 

S = \ J d 4 x v / — ^[Tfpji? - 2A] + S M , (1.3) 

where d 4 Xy/—g is the invariant volume element, g the determinant of the met¬ 
ric g^ u , R is the Ricci scalar and Mp\ is the reduced Planck mass with Mpj = 
(87TG)- 1 . 

By using the Principle of Least Action, this yields the Einstein field equations 

Gfiu + A = 8 TtGTfjv , (1-4) 


where G^ v is the Einstein tensor: 


G — Rfiv 2 Fadin' 


(1.5) 


It is composed of the metric g^ u , Ricci tensor R^ u and the Ricci scalar R = 
R a a = g al3 R a p. The right hand of this equation contains the matter within 
the spacetime, represented by the matter stress-energy-momentum tensor T^ u . 
The stress-energy momentum tensor is comprised of the energy density p and 
pressure p by the matter in question. It is defined as T ,l v = diag(— p,p, p, p) 


Eq. (1.4) describes the effect matter has on the curvature of the spacetime. 


In this Thesis, we will not be considering the Einstein field equations with 
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the inclusion of the cosmological constant A. 


A perfect fluid description of the system can be made. For this, the stress- 
energy momentum tensor can be written as 

T ^ = {p + p)u,j,u v + pg^ , ( 1 . 6 ) 

where u M is the 4-velocity of a comoving observer where 

dr} 1 

= — = 5^ , u^uT = -1 , (1.7) 

where A is an affine parameter associated with a geodesic in g^ v . 


Considering the FRW spacetime in Eq. (1.2) in general terms 


ds" = — d t 2 + a 2 (t)Sijdx l dx J , 


( 1 . 8 ) 


and applying this to the Einstein field equations in Eq. (1.4) and the energy- 


momentum tensor in Eq. (1.6), we obtain the Friedmann and acceleration equa¬ 
tions respectively 


3 H 2 = 8irGp , (1.9a) 

— 2H = 8nG(p + p) , (1.9b) 


where H = a/a is the Hubble parameter. Note that dots represents derivatives 
with respect to cosmic time. 
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The Einstein tensor fulfils the Bianchi identities and this results in V= 0. 
By taking the covariant derivative of the Einstein field equations, this leads to 
the conservation equations of the matter fields 


V /t T^ = 0 . 


( 1 . 10 ) 


For a system consisting of more than one fluid, the energy-momentum tensor will 
be composed of the different components; the total energy density is p = Yh Pi 

i 

and the total pressure is p = YPi- By taking into account an expanding 

i 

Universe (FRW spacetime) a conservation equation for each of the fluids can 
be derived 

Pi + 3H(l + Wi)pi = 0 , (1.11) 


where the parameter wy is the equation of state and is used to denote the matter 
in the system. It is given by 

w = P -, (1.12) 

P 

with some examples of the different equations of state shown below 


w = 


( 






0 

1/3 

-1 


matter-dominated , 
radiation-dominated , 
dark energy. 


We shall later study an extension of General Relativity in which the conservation 


equation in Eq. (1.11) is modified; in this model, there is an interaction between 


the scalar field and the matter species. 
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1.4 Inflation 

In the late 1970s, cosmologists realised there were problems with cosmology in 
its current state. They were encapsulated into three problems, these are the 
flatness, horizon and magnetic monopole problems. In Guth’s 1981 paper ra. 
the flatness and horizon problems were stated along with a mechanism that 
solves all these difficulties — inflation. 


1.4.1 The problems 

Flatness problem 

The flatness problem was stated by Dicke and Peebles in 1979 [60]. Through 
the observations made of the CMB, it is known that the energy density of the 
Universe is near critical density p cr i t . When the Universe is at critical density 
it means that the Universe is spatially flat i.e. k = 0. This can be shown by 
rewriting the Friedmann equation with the reinstatement of the curvature k, in 
terms of the density parameter, hi 


n - 1 


hi 


(1.13) 


with the parameter U defined as 



Pcrit 


(1.14) 
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and the critical density p cri t is given by 


Pcrit 


3 H 2 
8nG 


(1.15) 


The Universe has undergone both the radiation and the matter-dominated 
epochs, we shall consider both. By using the conservation equation stated 


in Eq. (1.11) and the relevant equation of state, a relationship between the 


energy density and scale factor is established, which is used in the first of the 
Friedmann equations to relate the scale factor to time. 


• Radiation-dominated epoch 

Using the method outlined in the section above, we can find relationships 
between energy density and scale factor and hence, the scale factor and 
time 

p ~ gU 4 =r- a ~ f 5 , (1-16) 


and use this in the rewritten Friedmann equation in Eq. 
radiation-dominated epoch where the equation of state is 
find 


(1.13). For a 


w 7 = 1/3, we 


aH 


~ t 2 , 


(1.17) 


and inserting this into the Friedmann equation yields 


— 1| ~ t . 


(1.18) 


• Matter-dominated epoch 


Following the same method applied to the radiation-dominated epoch, the 
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relationships between energy density and scale factor, and therefore the 
scale factor and time in a matter-dominated period with w m = 0 are 

p ~ a 3 ==>• a ~ f 5 , (1.19) 

and the resulting relation is formed 

aH ~ H , (1.20) 

with the resulting Friedmann equation 

\n - 1| - ■£§. (i.2i) 


In both the radiation and matter-dominated epochs, it is clear that |fl — 1| is 
an increasing function of time. From the results in the Planck 2013 paper EH. 
the total present day density parameter fi tot = ^ where i denotes the 

i 

matter fields such as radiation, baryons and dark matter, is 

fltot = 1-0005 ± 0.0033 . (1.22) 

This implies that 10 — 11 must be extremely close to zero in the early Universe, 
indicating that the initial conditions of our Universe must be incredibly fine- 
tuned. ffowever, such conditions are extremely unlikely, and cosmologists are 
continuously searching for dynamical reasons. 
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Horizon problem 


The horizon problem was first posed by Misner in 1967 [61]. We know from 
observations made of the CMB that the background radiation and therefore 


temperature, is approximately uniform in the Universe. From this, we can 
infer that regions on opposite sides of the Universe must have been in causal 
contact in the past. However, when extrapolating the light rays backward 
towards the surface of last-scattering, it is shown that these regions cannot 
have communicated. The surface of last-scattering is the point in time where 
the photons decoupled from charged matter. We now introduce the particle 
horizon (or comoving horizon) rj, and the comoving Hubble radius ( aH ) -1 , the 
particle horizon is defined as 



(1.23) 


and is the distance that light has travelled since the start of the Big Bang; 
a paper on the different horizons used in cosmology was presented in 1956 
by Rindler [62]. The comoving Hubble radius is the distance that determines 
whether two regions in the Universe can communicate or not. If we consider 
two particles that are separated by a distance greater than the horizon, there 
is no possibility of these particles interacting. However, if the distance is now 
greater than the comoving Hubble radius (and less than the horizon), the two 
particles could have possibly communicated in the past. The horizon problem 
involves answering why the CMB has a roughly uniform temperature, when the 
Universe contains many causally disconnected regions, approximately ~ 10 83 


regions. 
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Monopoles problem 

The problem of the magnetic monopoles was documented in 1978 by Zeldovich 
and Khlopov [63], and later in 1979 by Preskill [63]. It was restated by Guth 
and Tye [65] in 1979. In all the papers stated above, the Grand Unified Theory 
(GUT) was under consideration. The GUT is a theory that treats the electro¬ 
magnetic, string and weak forces as one force at very high temperatures, similar 
to those during the very early stages of the Universe. The aim is to combine 
this unified force with classical gravity in the hopes to understand the physics 
that takes place below the Planck scale. A feature of GUTs is that they predict 
the existence of stable magnetic monopoles at very high temperatures. 


1.4.2 The solution 


As a solution to three problems in cosmology, an inflationary period is required, 
which can be summarised in the relation 


(1.24) 


a > 0 . 


Using the definition of the Hubble parameter, the condition for inflation above 
can be rewritten in terms of the comoving Hubble radius: 



(1.25) 


Inflation is able to solve the flatness problem as it causes the term | O — 11 to drive 
towards zero, and hence present day density parameter stays approximately 
at unity. The horizon problem is solved as the comoving Hubble radius will 
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always decrease during an inflationary period, which allows more regions of 
the Universe to be in casual contact, and hence the CMB will be of uniform 
temperature. Due to the increase in volume generated by inflation, this causes 
any possible magnetic monopoles to be undetectable. 


1.4.3 Scalar fields and slow-roll inflation 

Let us now examine the second Friedmann equation. We can use this equation 
to study the accelerated expansion of the Universe; it can be written as the 
following 

- = -hl + 3w)p. (1.26) 

a b 

In order for this type of expansion to occur, a > 0, which can be obtained when 
the equation of state is w < —1/3. The scalar field is an example where this is 
possible, as we shall now show. 


The action for a canonical scalar field <fi with potential V (0) in Einstein gravity 
is given by 


S A 



Ig^d^d^-V^) , 


(1.27) 


and by varying this action with respect to the field, we can find its background 
equation of motion. The background equation of motion for the field <f> can be 
calculated by using the Euler-Lagrange equation 


SC f SC \ 


0 . 


(1.28) 
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Strictly speaking, a scalar field is dependent on both time and space i.e. 0 — 
x). However, we know that the Universe is both homogeneous and isotropic 
on large scales and therefore, we can say that the scalar field is dependent only 
on time at the background level. Using this statement and the Euler-Lagrange 


equation in Eq. (1.28) yields the Klcin-Gordon equation in FRW spacetime 


0 + 3 Hj> + U'(0) = 0 , 


(1.29) 


with primes denoting derivatives with respect to the associated field. 


The energy-momentum tensor of the scalar field is obtained by varying the 
action with respect to the metric g lllt 


2 SS 


= <9 m 0 d v (f) - g lw 


\g ap d a (f)dp(t) + V(4 >) 


(1.30a) 

(1.30b) 


By considering the leading diagonal of this tensor T ,J V = g^ a T au: the energy 
density and pressure of the scalar field are 


/+ = ^0 2 + V (0) , 

(1.31a) 

P* = ^0 2 - V (0) , 

(1.31b) 


and hence the equation of state for the scalar field is 

+ - 

w * if- + v(f ' 


(1.32) 
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The scalar field 0 is able to achieve the requirement for inflation, w < —1/3, by 
simply through the condition 0 2 -C V((f)). This is known as slow-roll inflation. 


Slow-roll inflation occurs when a scalar field (dominant for the case with more 
than one scalar field present) rolls down a potential U(0) [66]. In order to 
obtain this mechanism, the potential of the field must be significantly greater 
than its kinetic energy 0 2 -C V (0) and very slow acceleration of the scalar field 
is required 0 <C 1. Under these slow-roll conditions, the Friedmann equation 
and the Klein-Gordon equation for the scalar field are 

3 Hj> ~ -U'(0) , (1.33a) 

3 H 2 ~ U(0) . (1.33b) 


Adding to this, slow-roll inflation is usually defined using two parameters e and 
V 


e 


V 


H 2 2\Vj ’ 

0 V" 

#0 ~ Y ' 


(1.34) 

(1.35) 


These slow-roll approximations hold when e 1 and \r)\ 1. The formula¬ 

tion of the parameters were constructed by Licldle and Lyth in the early 1990s 
[671 68] . A formal approach of the slow-roll parameters is presented in [ 69] . 


A quantity that we shall require in the Thesis is the conformal time rj. The 
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relationship between conformal time and cosmic time is 



(1.36) 


During inflation, this unit of time will run from large negative numbers to zero. 
Another time variable which is useful in cosmology is the e-fold number N given 


by 



(1.37) 


where a enc j is the value of the scale factor at the end of inflation. 


1.5 Cosmological perturbation theory 

We will now move to the subject of cosmological perturbations. This area of 
cosmology has been studied extensively; see the following references for more 
detail m EH E3 [731GU- The linear perturbations of the metric g^ u , i.e. 
can be decomposed and classified depending on how they transform on spatial 
hypersurfaces; there are three types of metric perturbation: scalar, vector and 
tensor perturbations USES]. In this Thesis, we shall be studying two types of 
perturbations: scalar and tensor perturbations. 

Scalar perturbations are necessary in order to obtain structures in the Universe; 
some examples of structures include stars and galaxies. During the inflationary 
epoch, there were quantum field fluctuations that evolved to yield small density 
inhomogeneities in the Universe. This implies that there were seeds planted in 
the primitive Universe. The primordial density inhomogeneities were measured 
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to be approximately Sp/p ~ 10 -5 by COBE jl5]. Once the Universe reached 
the matter-dominated epoch, gravity causes these density inhomogeneities to 
grow in size leading to the formation of large scale structure. The other type of 
perturbation, namely tensor perturbations (or primordial gravitational waves) 
are a prediction of inflationary models, which can be used to eliminate such 
models. 

1.5.1 SVT decomposition 

The metric perturbations Sg /W can be decomposed into 4 scalar perturbations, 
2 vector perturbations and 1 tensor perturbation. We shall briefly summarise 
the scalar-vector-tensor (SVT) decomposition of cosmological perturbations. 

Scalar perturbations 

The line element for the scalar metric perturbations in the FRW spacetime is 

ds 2 = —(l+2'L)dt 2 + 2a(t)B i dtdx l +a 2 (t)[(l — 2$)5 i j + 2E i j]dx' l dx : ’ , (1.38) 

where T, B, $ and E are scalar perturbations of the metric. The gauge trans¬ 
formations can be made 


t —y t T ol , 

i' 4 i' + 8 lj /3j . 


where a is a scalar function and /3, is a divergence-free 3-vector. 


(1.39a) 

(1.39b) 
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Under these transformations the 4 scalar metric perturbations transform via 
the following 


T ->• 'I' - a , (1.40a) 

B —> B + a~ 1 a — a/3 , (1.40b) 

E ->• E-P , (1.40c) 

$ ->• + Ha . (1.40d) 


The matter perturbations given by Sp and Sp — which are the perturbations 
in the energy density and pressure — and the perturbation of the momentum 
potential Sq are also gauge-dependent and transform as 

Sp —> bp —pa , (1.41a) 

5p —>■ Sp — pa , (1.41b) 

Sq —» Sq + (p + p)a. (1.41c) 


Vector perturbations 

The vector metric perturbations S t and F t are given by this line element 

ds 2 = —d t 2 + 2aSidtdx l + a 2 [Sij + 2F^j)]dx l dx j , (1-42) 

where the perturbations satisfy the two conditions: S i:i = 0 and F l t = 0. The 
gauge transformation for these vector type perturbations are 

x l ->■ x l + P x , 


(1.43) 
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and vector type metric perturbations transforms as 

Si —> Si + afti , (1.44a) 

Fi^Fi-fa. (1.44b) 

Tensor perturbations 

The tensor metric perturbations h l3 are given by the line element 

ds 2 = —At 2 + a 2 [8ij + h i3 ]<lx l Ax 3 , (1-45) 

and hij can be written as follows 

K = h + ef 3 + h x e*j , (1.46) 

where ef 3 and ef 3 are the two polarization tensors and h + and h x are the 
coefficient functions. In addition, h i3 is both traceless and transverse US- 

5^hij = 0 , ffhij = 0 , (1.47) 

and the polarization tensors have the following properties 

e y (-k,A) = e*,.(k,A) and J]e*.(k,A)e«(k.A) = 4 , (1.48) 

A 

where A indicates the polarization i.e. A = + , x. Unlike the scalar and vector 
metric perturbations, the tensor metric perturbations are unaffected by gauge 


transformations. 
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1.5.2 Scalar modes 

In the previous subsection, we described the formalism required to study the 
first-order perturbations generated in cosmology. We shall begin with the scalar 
perturbations before turning our attention to tensor perturbations. There are 
two types of scalar perturbation that are produced during inflation; they are 
quantified as the curvature and isocurvature perturbations. 

Curvature (otherwise known as adiabatic) perturbations are those which are 
related to the perturbations in the total energy density (fptotai- The other type 
of perturbation is known as the isocurvature perturbation; this perturbation 
is one that does not affect the local curvature in such a way that it does not 
affect the total energy density perturbation 6p tota i = 0. For this perturbation 
type, the individual matter species components are allowed to vary. In the case 
of single-field inflation, there are no isocurvature perturbations generated. The 
presence of isocurvature perturbations causes growth in the curvature pertur¬ 
bation on superhorizon scales. 

In this Thesis, only the Newtonian (or longitudinal) gauge will be used. The 
name arises due to the system considered reducing to Newtonian gravity when 
in the small-scale limit; when at small-scales, the Newtonian potential T is that 
satisfying the Poisson equation V 2 T = 4nGp for nonrelativistic matter. In the 
Newtonian gauge, a coordinate transformation is chosen in such a way that the 
scalar metric perturbations B and E are equal to zero B = E = 0. This results 


in the full scalar perturbation FRW line element in Eq. (1.38) simplifying to 
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contain only $ and d/ which are now gauge invariant potentials 

ds 2 = -(1 + 2T)d t 2 + a?{t)( 1 - 2<f>)8 ij dx i dx j . (1.49) 

The perturbed Einstein field equations can be calculated by knowing the per¬ 
turbed Einstein tensors which are as follows 

5G% = — 2V 2< 3> + 6 H$ + 6 H 2 ^ , (1.50a) 

5G°i = -2(6 + Htyj , (1.50b) 

8G i j = [2$ + V 2 (T - <h) + (4 H + 6 H 2 )H + H( 2T + 6<h)]5} 

+ ^($-T). (1.50c) 

As an example, we can calculate the perturbed components of the energy- 
momentum tensor for a scalar field. We will assume that the scalar field is 
composed of two parts: the background (homogeneous) and the linear first- 
order perturbation 

x) = 0(t) + x) , (1.51) 

and we will be working with the Fourier components of the perturbations, Scj )k 

/ d 3 k 

(^^ kX > (i- 52 ) 

where Scp^ satisfies the Poisson equation: \ /2 5(pi i = —k 2 5(p. In order to shorten 
the expressions, the subscript k will be omitted. Furthermore, it is convenient 
to work with the gauge-invariant Sasaki-Mukhanov variables nznzgi. They are 
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related to the field perturbations and are defined as 


Q<t> = 50 + 



(1.53) 


From the perfect fluid description in Eq. (1.6) and the energy-momentum tensor 


for the field 0 in Eq. (1.30b), the perturbed components are 


n0 0 = -Sp = 0 2 ^ - 000 - V'8(j> , 

(1.54a) 

£T° =5q = , 

(1.54b) 

ST) =Sp= [-0 2 ^ + 050 + V'8<f>] 8) , 

(1.54c) 


where we have used 


= S(g^ a T au ) , 

= 6g^T au + g^6T a „ . 


(1.55a) 

(1.55b) 


In order to study the curvature and isocurvature perturbations, it is preferable 
to work with quantities that are gauge-invariant. For both types of perturba¬ 
tions in question a definition does exist for each. Starting with the curvature 
perturbations, first is the curvature perturbation IZ as defined: 

IZ = <F- 5q , (1.56) 

P + P 

8q is associated with the Oi-component of the perturbed Einstein field equations. 
We can show that this curvature perturbation is a gauge-invariant quantity. By 


using the definition in Eq. (1.56) and considering the scalar metric perturbation 
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gauge transformations in Eqs. (1.40) and (1.41) 


IZ — ^ + Ha - \5q+(p + p)a] , 

p + p l ' J 

= $ + Ha - 5q — Ha , 

P + P 

„ ^ H 

1Z = - 5q . 

P + P 


This type of curvature perturbation represents the gravitational potential on 
comoving hypersurfaces where the inflation fluctuations are 5<p = 0 


n = 4>| 50=o • (1-57) 

For the case of a single scalar field, the comoving curvature perturbation is 

K = $+S-6<l>. (1.58) 

<t> 

Second is the curvature perturbation on uniform density hypersurfaces £ m 

- C = $ + — $p , (1-59) 

P 


and is gauge-invariant (using Eqs. (1.40) and (1.41)) as shown below 









32 


CHAPTER 1. INTRODUCTION 


This curvature perturbation describes the gravitational potential on uniform 
energy density slices: 

- C = ®\s P =o • (1.60) 


Using the conservation equation in Eq. (1.11), the expression for ( can be rewrit¬ 
ten as 


-c = $- 


5p 


3 (p + p) 


(1.61) 


During the slow-roll period in the case of a single scalar field, the following hold 
true 


p + p = (f> , 

5p ~ —3 H<fi8(j) . 


(1.62) 

(1.63) 


The last relation conies from using the Klcin-Gordon equation of the scalar 
field; first the field perturbation 8<p is frozen in and under slow-roll conditions 


Eq. (1.33a) is obeyed. This results in the following relation at superhorizon 


scales 


-c = * 


3 Hcj) 
3 0 2 


_ H 

= + —8(j) , 


(1.64a) 

(1.64b) 


which is the same as Eq. (1.58). Therefore at superhorizon scales, both — £ and 
1Z are the same. 


The other type of perturbation that can be calculated is the isocurvature per- 
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turbation. For a system to be adiabatic, the pressure perturbation is as follows 

dp = c^dp , (1.65) 

where c ^ the adiabatic sound speed of the fluid given by = p/p. However, in 
general, the total pressure perturbation is expressed containing two contribu¬ 
tions: an adiabatic part and an entropic part also known as the nonadiabatic 
pressure perturbation dp nad m- 

Sp = clip + <5p„ad , (1.66) 


where <5p na d is defined as dp na d = pT. The entropy perturbation T is defined as 


v = 6 P_^P 
P P 


(1.67) 


which is by construction a gauge-invariant quantity as shown below 


5p — pa 5p — pa. 

p p 

dp dp 

— — — a — — T a , 
P P 

p _ 5 JP_ _ 

P P 


This quantity describes the difference between hypersurfaces of uniform pres¬ 
sure and uniform density [8Uj. 


An isocurvature perturbation that is regularly calculated in literature is known 
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as the entropy perturbation S which is defined as [8IJ E2J 

S = — hp nad . (1.68) 

P 

Earlier studies into the entropy and nonadiabatic perturbations are shown in 
references [[83] and [M] . The nonadiabatic pressure perturbation also appears 
in the rate of change of the curvature perturbation on hypersurfaces of uniform 
density; the full expression is presented in [79] 

C =-flpnad - \v 2 (o + v + B) , (1.69) 

p + p 3 

where VT is the perturbed 3-velocity of the fluid, a denotes the shear given by 

a = E — B , (1.70) 


where B and E are the metric perturbations in the line element stated in 


Eq. (1.38). At large scales, the gradient term in Eq. (1.69) can be neglected 


resulting in 


c = 


H 


P + P 


$Pna,& ■ 


(1.71) 


We require the primordial power spectrum as it can be used to distinguish 
between various inflationary models. Some of these models produce features in 
the primordial power spectrum in the form of bumps and oscillations, examples 
of models include those with features in the effective inflaton potential [ 851 
m E3 m Ha [90] and disruption to the slow-roll evolution caused by phase 
transitions [91]: for more examples see [92]. 
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The power spectra for the curvature and entropy perturbations and their cross 
spectrum between the two quantities are stated below 

k 3 

V ( = — 2 m 2 + \C 2 \ 2 ), (1.72a) 

k 3 

V s = — 2 (\S 1 \ 2 +\S 2 { 2 ) , (1.72b) 

k 3 

Pc = Ci £i + C 2 $2 I > (1.72c) 

where k is the wavenumber, this quantity arises in the perturbation equations, 
and the subscripts indicate the run number. For an overview of the numerical 
method performed in this Thesis, see Appendix A. 


Another quantity that is commonly seen with the cross spectrum is the cor¬ 
relation rc which is defined as 


Pc 

Tr — - 

s/Ws 


(1.73) 


and varies between the values zero and one. Once the power spectrum has been 
obtained, other quantities can be derived from it. The power-law spectrum of 
the primordial curvature perturbations is defined as 


*><(*) = w(£) 


( 1 . 74 ) 
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where V^(k 0 ) is the amplitude of the curvature power spectrum at a pivot scale 
denoted by k 0 and the spectral index n s is 


n s - 1 = 


d In Vq 
d In k 


(1.75) 


Some models predict a scale dependence, this is also known as the running in 
the almost power-law spectrum of curvature perturbations. For this case where 
the power spectrum contains running [93], it is given as 

( , \ n s —1+| ln(k/k 0 )a 

h) ’ ( L76 ) 

with a as the running of the spectral index. 


dn s 
din A: 


(1.77) 


In the WMAP experiment, the pivot scale is set at ko = 0.002 Mpc 1 |9 ij . 


1.5.3 Fluctuation origins 

It was proposed in early 1980s by various theorists that scalar fields are sub¬ 
jected to quantum fluctuations with these fluctuations stretched when on su¬ 
perhorizon scales [in, 1H EH EM; this was also looked at later in [1U0] . The 
quantum field fluctuations, which originate at very small scales, will evolve dur¬ 
ing inflation and will become the initial density perturbations which will give 
rise to structures we see today. 
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During the inflationary period, the fluctuations of the inflaton field <f> will grow 
until a point where its wavelength is greater than the Hubble radius. The am¬ 
plitude of the field fluctuations will then freeze at a nonzero value and will 
remain so until the end of inflation. At some stage, inflation eventually ends 
causing the field fluctuations to reenter the Hubble radius due to the Hubble 
radius increasing faster than the scale factor. The reentering of the fluctuations 
does not happen immediately, but occurs after reheating during either the ra¬ 
diation or matter-dominated epochs. For the fluctuations that exit the horizon 
between 50 and 60 e-foldings before the end of inflation, these will reenter at a 
scale corresponding to cosmologically relevant wavelengths. When these modes 
reenter the Hubble radius, they will leave an imprint in the curvature spectrum 
which can be measured and compared to various models. 


The fluctuations of the scalar field can be studied by using its second-order 
action; the second-order action of the scalar field (j) written in terms of the 
canonically-normalised Mukhanov variable u in conformal time 77 is given by 


S® 



( diU ) 2 + —u 2 

z 


(1.78) 


where the parameters u and z are as follows 
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where IZ is the comoving curvature perturbation. This calculation is shown ex¬ 


plicitly in the following references: mmm- Varying the action in Eq. ( |1.78[ ) 
yields the Sasaki-Mukhanov equation 


“k + w k (>j) “k = o , w k = ( k - — ) , 


(1.80) 


with the following Fourier modes defined as 


«k(>)) = / d 3 xe ‘ kx a()),x) . 


(1.81) 


In a de Sitter spacetime which is a spatially flat spacetime containing a positive 
cosmological constant, the scale factor is 


a{v) = 


Hr] ’ 


(1.82) 


which yields an effective oscillator frequency of 


2 i 2 * 

uj \, = k - — 2 


(1.83) 


This example in de Sitter spacetime can also be studied in two different regimes: 


• Subhorizon scales 

Subhorizon scales occur when the wavelengths of the fluctuations are much 
smaller than the Hubble radius A H and hence k 2 \z"/z\, leading 
to a frequency of 


2 / 2 
ui k = k , 


(1.84) 
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and yielding the oscillatory solutions to the Sasaki-Mukhanov equation 


°c e 


zLikrj 


(1.85) 


• Superhorizon scales 

Superhorizon scales are when the wavelengths are much greater than the 
Hubble radius A H 1 and therefore k' 2 <C \z"jz\ 



The resulting solution for this regime is 


( 1 . 86 ) 


oc - and n k oc 77 
V 


(1.87) 


Using the fact that during inflation conformal time will run from large negative 
values to zero, we have here the first of these solutions is a growing mode and 


the second is a decaying mode. From using Eq. (1.79a) and that z ~ 77 , it 

can be shown that the comoving curvature perturbation 77 is constant on su¬ 
perhorizon scales. 


A vacuum state will need to be chosen and the standard choice is to use the 
Bunch-Davies vacuum, named after the authors whose work was published in 
1978 [Ml. At very early times, the wavelengths of all fluctuations existed in¬ 
side the Hubble radius (subhorizon scales) and therefore the following equation 
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holds 

«k + k 2 u k = 0 , 


( 1 . 88 ) 


and through the quantization of the field u k the solution to the equation above 

is 

lim u k {rj) = — ^=e~ lkri . (1.89) 

v^-oo v2 k 

The power spectrum for the scalar modes is 


Vn 


k 3 

Mv) 2 

2 t r 2 

z(v) 


k 3 
2 ^ 


l^ k (h)| 2 ■ 


(1.90a) 

(1.90b) 


1.5.4 Tensor modes 


Inflation also predicts the existence of gravitational waves which are tensor 


modes. The line element for tensor perturbations is given in Eq. (1.45) We 


would like to study the fluctuations of the tensor models and in order to do 
this, we will need to find the second-order action of the tensor modes. This 
is done by expanding the gravitational sector of the action (Einstein-Hilbert 
action) to second-order; the second-order action for the tensor perturbations is 


= I dr?d 3 x^ 


{h'Y - (dihy 


(1.91) 


which will lead to the equation of motion for the quantity /ij. 


/?k + 2 — h k + k 2 h k = 0 . 


(1.92) 
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This equation can be rewritten by defining a new variable Uk where 

t>k = ahk , (1.93) 


thus yielding the equation 


a 

a 


v k + ( k - — ) v k = 0 . 


(1.94) 


It is important to note that this equation is of the same form as the Sasaki- 


Mukhanov equation for scalar modes in Eq. (1.80), with the exception of the 
scale factor replacing the parameter z. 


The same arguments used on the scalar modes can be applied to tensor modes 
- all wavelengths in the very early Universe existed within subhorizon scales. 


Due to this, the term a"/a can be neglected and Eq. (1.94) is reduced to that 


in Eq. (|1.88|). The initial conditions for the gravitational waves are the same as 

(1.95) 


for the scalar modes 


t’k = 


D —ikri 


V2k 


The power spectrum generated for the tensor modes is 


V h = 


k 3 
2 ^ 


Mv) 


a (v) 


(1.96a) 


= ■ 


(1.96b) 









42 


CHAPTER 1. INTRODUCTION 


In order to relate this to the gravitational waves (tensor) power spectrum, we 
will need to sum over all polarizations (appears as the additional factor of 8): 


Vt 


4x2 


k 3 

2^ 


\hM \ 2 , 


= sv h , 


(1.97) 


A quantity that is commonly used to relate the tensor power spectrum Vt to 
the scalar curvature power spectrum V\ is the tensor-to-scalar ratio, r and is 
defined as 


Vt _ 8 Vh 

vl~l\ 


(1.98) 


1.6 Alternative models 


Most inflationary models are based upon scalar fields, both single and multifield, 
in General Relativity. However, modifications to Einstein’s Theory of General 
Relativity have also been studied, for example through scalar-tensor theories, 
such as Brans-Dicke theory |102] and f(R) theories |T03| Il04j . For further 
reading on f(R) theories see the following reviews [ 1051 196] . and a more general 
overview of the various extended models of gravity see the review m 
The action of a scalar-tensor theory is usually displayed as 


S= d 4 x^g 


^F^R-^d^d^-UicP) 


+ Sxlg^x] > (1.99) 


where F(0) is the coupling between the scalar field and the gravity sector. The 
action S x [g IJLI ; x] is that of the matter fields y; this action is only composed of 
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the metric g^ v and the field x and is independent of the field 0. By rewriting the 
action, the Einstein-Hilbert term can be recovered. The most common practice 
is to perform a conformal transformation of the metric given by 


9/j.i/ R i.4*) 9 , 


( 1 . 100 ) 


where g lw is the redefined metric in the new representation, and a scalar field 
redefinition detailed below 


di>y 3/F 0 y i 
d(t>) 2\F ) F ' 

(1.101) 

a($) = f-*{<i>) , 

(1.102) 

V($) = f/(0)^ 2 (0) . 

(1.103) 


By using the new definitions stated above, the action in Eq. (1.99) is rewritten 
as follows 


S = / dW- 




+ S x [A 2 (4>)g^;x] i (1.104) 


where g is the determinant of the metric g lw and R is its Ricci scalar obtained 
by using the tilde metric. The usage of the conformal transformation in scalar- 
tensor theories are presented in the following references: [108 .. 1109] . 


In the presence of a second scalar field, these theories will provide the sys¬ 
tem with nonstandard kinetic terms, such as kinetic couplings between scalar 
fields. 
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The idea of a Lagrangian containing nonstandard kinetic terms coined as “k- 
inflation” was first introduced in the papers [110 ' . Illlj published in 1999 as a 
possible solution to inflation, using inspiration from models in string theory. 
From string theory based models, an inflaton candidate known as a moduli 
field which is a weakly coupled scalar field, naturally arises. However, it was 
shown that the potentials generated from these models do not allow for slow 
roll to occur. It was later shown that these nonstandard kinetic terms are able 
to produce interesting models in the context of dark energy; the class of models 
called “/c-essence” mat mm-, for an extensive review of dark energy see 

mg. 


The kinetic couplings between scalar fields are obtained by performing the con¬ 


formal transformation on the action stated in Eq. (1.99); works in this area 


include mg nn ng mi dan . 


1.7 Thesis outline 

This Thesis is divided up into five chapters, with the first chapter dedicated to 
the history of cosmology and the mathematical formalism required. The follow¬ 
ing three chapters will each discuss models involving coupled scalar fields under 
a different setting; Chapters 2 and 3 will consider the inflationary scenario and 
Chapter 4 in the context of dark energy. 

Chapter 2 considers three models of multifield inflation each differing in ki¬ 
netic terms, and their effects upon the isocurvature (nonadiabatic pressure) 
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perturbations generated during the inflationary epoch. Each of the models con¬ 
sidered contains at least one minimally conpled scalar field, with two of these 
models classified under /e-inflation. The models studied are as follows: the first 
model containing two minimally coupled scalar fields, second model consisting 
of two fields, one minimally coupled and the other containing a kinetic cou¬ 
pling, and lastly, the third model composing of a minimally coupled field and a 
Dirac-Born-Infcld field. 

In Chapter 3, a study of the effects induced by sharp transitions in the effective 
Planck mass upon primordial power spectra during inflation is presented. A 
single field model of scalar-tensor form is initially studied before extending to 
a two-field model with the inclusion of a minimally coupled auxiliary field. In 
both cases, the scalar and tensor power spectra are generated numerically. It 
was shown that the sharp transitions affected slow-roll evolution, and this was 
reflected through a feature in the scalar power spectrum. This feature allowed 
the easing of tension between observations from the two experiments, Planck 
[57] and BICEP2 [121] . 

The move to more later times in cosmic history occurs in Chapter 4 where 
a noncanonical transformation is studied in the context of dark energy. In this 
chapter, a model of scalar-tensor theory form is presented in which the scalar 
and matter fields are disformally coupled. The intention is to study the effects 
disformal coupling has upon observables, such as the temperature and spectral 
distortions of the CMB. In order to check the consistency of the disformal theory, 
a fluid description and a kinetic theory was verified. Through the verification, 
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two important traits of the disformal theory were discovered: the equation of 
state and the distribution function are not frame-independent quantities. 

The Thesis concludes with a final chapter recalling and summarising the main 
points from the three studies contained in the preceding chapters. Suggestions 
to extend the work performed in this Thesis and concluding statements are also 
included. 

Notation 

We will not be considering the cosmological constant A. For the most part of 
the Thesis, the reduced Planck mass Mpi will be set to one Mpi = 1, except in 
fundamental equations such as the action of a theory. All the calculations in 
this Thesis will be performed in the Newtonian gauge. 

The metric signature used throughout this Thesis is (—,+,+,+). For mathe¬ 
matical objects such as the metric g^, Greek indices i.e. (a, (3 ,...) run from 0 
to 3, and denote over all dimensions, and lower case Roman indices i.e. (a, b ,...) 
run from 1 to 3, and represent the spatial coordinates. 



Chapter 2 


Noncanonical multifield inflation 


The principles of inflation were laid out in the Introduction, alongside a brief 
note of alternative models in its self-titled section. There are numerous theories 
beyond the standard model of physics that provide motivation for inflation¬ 
ary models; some examples of these theories include those from supergravity 
[ 1221 123 . 1124] and string theory [12511126] . As a result, these models predict 
many scalar fields that could drive the inflationary epoch. In addition, these 
scalar fields couple through either the potential or kinetic terms. We shall be 
concentrating on the details of the latter. 

Many of these extension of standard physics models give rise to fields with non¬ 
canonical kinetic terms; this group of models are known as 1-inflation. Some 
examples of work in this field include [1101111111127111281112911180111811118211188] . 
One example of 1-inflation is one driven by the Dirac-Born-Infeld (DBI) field, 
which has been studied in [4341 435] . 
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Inflation in the DBI model is caused by a D3-brane travelling along a warped 
region called a throat of a compactified space. Additionally, a speed limit 
naturally arises on the brane restricting the speed at which it travels, and is 
dependent on the speed of the brane and the warping of the throat. With this 
speed limit, a parameter is introduced which is the analogue of the Lorentz fac¬ 
tor in special relativity and is allowed to grow until the speed limit is reached. 
This boost factor also has an effect on the perturbations generated. Unlike 
fluctuations of standard scalar fields which travel at the speed of light, the 
fluctuations in the DBI field travel at a sound speed which is related to the 
analogous Lorentz factor. 


Features from these noncanonical models with multiple fields can be used to 
distinguish and eliminate models that do not agree with current constraints 
set by observations of the Cosmic Microwave Background (CMB) radiation 
[ 136 1 . T3? l I138j . One such constraint is through using isocurvature perturba¬ 
tions; the amount of isocurvature perturbations is restricted to approximately 
of order 10% of the curvature perturbations [941j . 


Recently, the nonadiabatic pressure perturbation has been identified as an im¬ 
portant ingredient in the study of cosmological perturbation theory. The per¬ 
turbation can source vorticity perturbations at second order in cosmological 
perturbation theory j 139j . In fact, the second order vorticity is shown to evolve 
as 


u 2ij ~ 3Uc 2 s uj 2ij 


2 a 

P + P 


3F^Upi^Pnadl ,j] + 


& Pl,\j&Pnaxil,i] 


P+P 


( 2 . 1 ) 
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where 5p\ and <5p na di are the first order energy density and nonadiabatic pres¬ 
sure perturbations. It was found that if the first order nonadiabatic pressure 
perturbation vanishes, it causes the vorticity to decay at both first and second 
orders. In contrast, under the context of multi-field inflation, the nonadiabatic 
pressure perturbation is usually nonzero and therefore sources vorticity at sec¬ 
ond order m um eh . Vorticity can possibly be detected by searching for 
B-rnode polarization of the CMB radiation [142] . hence creating further tests 
of inflation. 

The evolution of the nonadiabatic pressure perturbation has been studied in 
detail, specifically for the theory involving two fields with canonical kinetic 
terms and a variety of potentials by Huston and Christopherson in [143] . The 
potential choices included the double quadratic m , quartic [1451 1146] and 
product exponential [22]. In their paper, the authors stated and used two 
known definitions of the entropy perturbation: one definition in which directly 
uses the nonadiabatic pressure perturbation, and the other more commonly 
used definition of utilising field decomposition as proposed by [51] . Both these 
entropy perturbation definitions were used for each potential considered. With 
each potential, the evolution of the power spectra with respect to the e-fold 
number and wavenumber for the curvature and entropy perturbations were 
plotted. From this, they found that the entropy perturbation which uses the 
nonadiabatic pressure perturbation evolves differently compared to that using 
the regularly used isocurvature mode through field decomposition. 


In this chapter, we extend the works by Huston and Christopherson in 2012 
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m , by studying various models with noncanonical kinetic terms, concen¬ 
trating on the method which incorporates the nonadiabatic pressure pertur¬ 
bation. In particular, we will consider two theories in which both will contain 
a canonical kinetic field: one theory with the inclusion of a kinetic coupling 
term between the two fields [1191. 1481. Il'4'9j , and the other with a DBI field, 
see [1341 11351 1T5UI :151] and references therein. This chapter is presented as 
follows: the actions for the three models considered are stated in the following 
section along with their background and perturbation equations in Section 2.1[ 
Section [212] presents the results for the three models considered. For each model 
apart from double inflation, there are two scenarios: one in which the canonical 
field dominates the inflationary epoch, and the other, vice versa. A summary 


of the findings are presented in Section 2.3 


2.1 The models 


We will consider three models: one model containing two canonical scalar fields, 
the second containing two scalar fields with a kinetic coupling between them 
and lastly, a model containing a canonical scalar field and a DBI field. Their 
actions are given as follows: 

1. The first two models are described by actions of the form 


S 



r - \<rd^dA 




( 2 . 2 ) 


with &(</>) denoting the kinetic coupling between the two fields (j) and x- 
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In the first model, there is no kinetic coupling and so b((f>) = 0, whereas 
in the second model b is defined as 

b(<i>) = M (2.3) 


where [3 is a constant. 


2. The third model we examine contains a scalar field with a canonical kinetic 
term 0, and one DBI field y. The action is given by 


S= d 4 x^/^g 


Mk R _ 




fix) 




, (2-4) 


where the two parameters exclusive to DBI field are defined by the fol¬ 
lowing 


7 = 

fix) = 


1 

V 1 + fix^dvXduX ’ 
A 

ix 2 + h 2 ) 2 ’ 


(2.5a) 

(2.5b) 


7 is a parameter which is analogous to the Lorentz contraction factor in 
special relativity and is related to the sound speed of the DBI field by 


( 2 . 6 ) 


It is important to note that c s is the speed at which the fluctuations of 
the DBI held y propagate, /(y) is the warp factor which describes the 
shape of the extra dimensions. The two parameters A and /j in the warp 
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factor are constants. 

In all the actions given, Mp\ is the reduced Planck mass where Mp t = (87 tG) _1 , 
R is the Ricci scalar and V(eft, y) is the potential. We shall set the reduced 
Planck units to one, Mpi = 1. 


2.1.1 Background 


We assume a spatially flat expanding spacetime which is both homogeneous and 


isotropic. This spacetime is given by the FRW line element as given in Eq. (1.8) 


in Section 1.3 There are no modifications to the Einstein-Hilbert part of the 
three actions and hence, all the actions obey the Einstein field equations 


— T^ v 


(2.7) 


where is the Einstein tensor and T tw is the energy-moment urn tensor. We 
assume that the energy-momentum tensors of the models considered are of a 
perfect fluid form, T^ u = diag(— p,p,p, p) where p and p are the energy density 
and pressure. 


The background Klein-Gordon equations are obtained by varying the actions 
with respect to the fields (ft and y. Due to no change in the form of the Einstein 
field equations, the Friedmann and acceleration equations are the same as those 


in Eq. (1.9) 


For the actions stated in Eqs. (2.2) and (2.4), the Klein-Gordon equation, energy 


density and pressure are summarised in the following: 
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• Kinetic coupling 

The Klein-Gordon equations are found by varying the action with respect 
to the two fields and are as follows: 




V ft (^TV^) = -b^V^uX ~ V tp , 




V M (^^e 2ft V, X ) = V, 


x • 


(2.8a) 

(2.8b) 


The FRW metric is applied to the equations above and yield the equations 
of motion in cosmic time. In addition, the energy density and pressure, 
and hence Friedmann equations are also stated here; they are given by 

maun] 


(f> + 3 Htj) + 14 — b<f,e 2b x 2 , (2.9a) 

X + (3 H + 2b ( j > 4>)x + e 2b V x = 0 . (2.9b) 

The energy-momentum tensor is given by 

T, 1V = - g, lv Q^ a/3 0, Q 0,/3 + v'j 

+ e 2b X,nX,v * 9^ (^9 af>e2b X, a X,pj • ( 2 - 10 ) 


Both the energy density and pressure are obtained by considering the 
leading diagonal of the energy-momentum tensor (the same method as 
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employed in Section 1.4.3) 



(2.11b) 


(2.11a) 


where b$ = db(4>)/d(j); subscripts denote differentiation with respect to the 
subscript. For all the equations stated above, they will reduce to the two 
standard field case when the kinetic coupling is set to zero. 

• DBI field 

The Klein-Gordon equations for the two fields in the DBI model are 


-j=y,{V^~gg^A) = Vj, , 


(2.12a) 



(2.12b) 


and when the FRW metric is applied it yields 


cp + 3H(f) + V ( p = 0 , (2.13a) 


X + 3^7“ 2 X + ~ 37" 2 + 2 7 " 3 ) + 7" 3 ^ = 0 , (2.13b) 


where f x = df jd\. The energy-momentum tensor for this model is 



(2.14) 
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From this, the energy density and pressure are 

P = + y( 7_ 1 ) + 1/ ’ (2.15a) 

p=\i 2 + l(l-L)- v . (2.15b) 

From the definitions above, the Friedmann equations are 

%H 2 — + )(7 — 1) + V , 

-2// = -r-':V , 

where the following relation 

1 ~~ 2 =fx 2 (2-17) 

7 

has been utilised. Likewise to the case containing the kinetic coupling, 
when the boost factor is set to one, the equations reduce to those of the 
two canonical fields. 


(2.16a) 

(2.16b) 


2.1.2 Perturbations 


We turn now our attention to the first order perturbation equations. First, we 
choose a gauge to perform this calculation in; for this work, we will perform 


the calculation in the longitudinal gauge, which was discussed in Section 1.5.2 


In this gauge and without the presence of anisotropic stress, the two scalar 
metric perturbations 'k and are equal. The perturbed FRW line element is 


essentially the same as that given in Eq. (1.49) but now including = T, we 
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will repeat it here for convenience 


ds 2 = - (1 - 2T)d t 2 + a 2 (l + 2'!/)8 ij dx i dx j . (2.18) 


The perturbed Einstein field equations in the longitudinal gauge are 


k ' 2 

3 H(HH + T) + —T = 
a 2 


(2.19a) 

T + TTT = 

~l Sq > 

(2.19b) 

T + 4TTT + (2W + 3TT 2 )T = 


(2.19c) 


where the perturbations in the energy density, momentum potential and pres¬ 
sure are dp, dq and 5p, respectively, as stated in Section [l.5.2| We will now state 
the perturbation equations — both the Klein-Gordon and Einstein equations 
— for the models containing the noncanonical kinetic terms. 


We begin with the model involving the kinetic coupling between the two fields. 
The components, dp, dq and dp, of the perturbed energy-momentum tensor are 


calculated using Eqs. (1.55) and (2.10), and are as follows 


dp = <j)d<j) + e 2b \d\ + b^e^x 2 dcj) + V^dcj) + V x d\ - - e 2b y 2 T , (2.20a) 

dq = 4>d(j) + e 2b xdx , (2.20b) 

dp = <j)d<j) + e 2b xdx + bcj,e 2b x 2 d(f> — V<pd(j) — V x dx — </> 2v k — e 2 b x 2 T . (2.20c) 
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The perturbation equations for the two fields are given by perturbing the general 


form of their equations of motion given in Eqs. (2.8), and are given as 


<50 + 3 77<50 + 


k 2 


'—^+V H - (b H + 2 &J)0 2 e 2b <50-2 ^e 26 y<5y 

CL 

+ V<Px S X ~ 4 # + 2V^ = 0 , 

$X + (3 H + 2 b^4>)8x + \ + e~ 2b V xx 5\ + 2^x<50 

CL 

+ e 26 (V^ — 2b ( f > V x ) + 2b ( j )( j ) 4>x 


(2.21a) 


<50 - 4\T + 2e~ ib V0'b = 0 . (2.21b) 


We can rewrite these perturbation equations in Eq. (2.21) in a gauge-invariant 


form using the Sasaki-Mukhanov variables stated earlier in Eq. (1.53) in Sec¬ 


tion 1.5.2 they are given as follows 
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with the coefficients C^, C<p x , C xx and C X( p as 
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The perturbation equations (Klcin-Gordon and Einstein equations) for the 
model involving the kinetic coupling between the two fields can be found in 
the papers [TT91 I12UI1152] . 


We now move to the perturbation equations for the model concerning the 
canonical scalar and DBI fields. In this model, the energy density, momen¬ 
tum potential and pressure perturbations are as follows 

5p = 050 - 0 2 T + 1050 + V x Sx + “ ^2 - 3y + 5\ 

-7 3 (X<5X + X 2 ^) , (2.24a) 

5q= - (050 + 7X<5x) , (2.24b) 

dp = 050 - 0 2 T - 10,50 - V x Sx 6x 

+ l(x$X ~ X 2 ^) • (2.24c) 


We find the Klcin-Gordon equations for the field perturbations are 


50 + 3#50 + ( % + T0J50 + VfrSx ~ 404 / + 2V^ = 0 , 
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By rewriting the field perturbation equations in Eq. (2.25) using the Sasaki- 


Mukhanov variables as stated in Eq. (1.53), the gauge-invariant form of the 
perturbation equations are 


Q(/> + 3 HQ'f, + B ( j ) Q x + + C ( f >x Q x — 0 , (2.26a) 

Qx + + 3-^ Q x + 5 X Q + C xx \ Qx + C x4>Q* = 0 5 (2.26b) 


with the coefficients BB x , C^, C<p x , C xx , C X 0 , in the equations are as follows 
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Further on in this work, we shall calculate the power spectra for both the curva¬ 
ture and isocurvature perturbations. We will now recall elements of the different 


perturbations from Section 1.5.2 


The curvature perturbation that we have chosen to study is the curvature per¬ 


turbation IZ in Eq. (1.56). It can be rewritten in a different form by using 


the acceleration equation in Eq. (1.9b) and the perturbed Einstein equation in 

H 2 ( T\ 

R=--g-U + -J+<I. (2.28) 


Eq. (El9bk 


For the isocurvature perturbation, we shall calculate the entropy perturbation 
S and this is a gauge-invariant quantity [HUES], an d if i s related to the nonadi- 
abatic pressure perturbation 5p na d. The gauge-invariant entropy perturbation 
and the nonadiabatic pressure perturbation are related through the definition: 


<S = —5p nad . (2.29) 

P 

The nonadiabatic pressure perturbation is a constituent of the pressure pertur¬ 
bation 5p] the pressure perturbation can be written as 


Sp = c 2 Jp + 5p nad , 

c2 = i 

s,a p l 

with c s 2 a denoting the adiabatic sound speed. 


(2.30) 

(2.31) 


The adiabatic sound speed for each of the three models can be calculated by 






2.1. THE MODELS 


61 


using the time derivatives of the Friedmann and acceleration equations, and the 
Klein-Gordon equations for the two fields. We now state the adiabatic sound 
speeds for the three models considered in this work. This speed should not be 
confused with the sound speed c s 2 of the DBI held fluctuation. 


There are two speeds that we use in this work: the adiabatic sound speed 
c s 2 a and the DBI held perturbation sound speed c 2 . To be clear, we now make 
the distinction between the two. The hrst of these two speeds, the adiabatic 


sound speed, is given by the standard definition stated in Eq. (2.31), and the 


second speed is introduced when considering the speed at which a held pertur¬ 
bation propagates. This is the held perturbation sound speed, otherwise known 
as the phase speed; for further discussion on this topic, see [153j. 

1. Two canonical scalar fields 


2 — 1 i 2 ( V <P^ + V xX) 

a. -L ~T 


3H(cj) 2 + x 2 


(2.32) 


2. One canonical scalar field and one scalar field with a kinetic coupling 


c 2 _ 1 + V X X) 

3 H^ 2 + e 2b x 2 ) 


(2.33) 


3. One canonical scalar field and one DBI field 
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For the final two models that contain noncanonical kinetic terms, the adiabatic 
sound speed reduces to the two canonical scalar field model when the kinetic 
coupling and boost factor are b = 0 and 7 = 1. 

Due to the complexity of these equations, for all models considered, we will 
solve them numerically following the method outlined in [120, T54i 1T5'5] . 


2.2 Results 


We will now describe the results of our numerical calculations. To be con¬ 


crete, we consider all models presented in Section [271] with the double quadratic 
potential [144] which will be rewritten in the form 


mx) = ^(f+ r V), 


(2.35) 


with T as the ratio between the two field masses 


m. 


p _ ■•"x 

m<p ’ 


(2.36) 


where rri^ and m x are the masses of the fields (j) and y, respectively. 


All the plots in this chapter are shown at the Wilkinson Microwave Anisotropy 
Probe (WMAP) pivot scale where &wmap = 0.002 Mpc^ 1 [94j and at plotted 
against the number of e-folds. For all models considered, the parameter values 
have been chosen so that the final amplitude of the curvature power spectrum 
is Vtz ~ 2 x 10 -9 at the WMAP pivot scale. 
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Figure 2.1: Background field evolution for the model containing two canonical 
scalar fields: 0 (solid red) and x (blue dashed). The initial field values are 
0 ini = Xi n i = 12.0 and masses of the two fields are V = 7.0 and = 1.395x 10 -6 . 


2.2.1 Two canonical scalar fields 


In this model, we consider the case where T = 7.0 as studied in [ 1201 HT5j . To 
further this work, we impose a further condition. This condition is to match 
the curvature power spectrum amplitude of approximately V-r. ~2x 10~ 9 at 
the WMAP pivot scale. In order to do this, the mass required by the field 0 is 
= 1.395x 10 -6 . The initial field values for the two fields are y ini = 0 ini = 12.0 
as in [143]. In addition, the two field derivatives at the start of inflation are 
those in slow-roll. The background dynamics for this model are presented in 
Figure |24~ and the comparison of the various power spectra evolutions, includ¬ 
ing the curvature, entropy, pressure and nonadiabatic pressure perturbations, 


are displayed in Figure 2.2 


We see in Figure |2.1 that the field y reaches the minimum of the potential, 
resulting in the 0-field dominating for the last 30 e-folds of inflation. This is 


reflected in Figure 2.2, which at this same point in time, there is a rise in the 
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Figure 2.2: Comparison of the various power spectra evolution with respect 
to e-fold number generated for the two canonical scalar field model. The initial 
field values are d> ini = Xmi = 12.0 and masses of the two fields are T = 7.0 and 
= 1.395 x 10~ 6 . 

Left panel: Power spectra evolution for the curvature (Vn, blue dashed) and 
entropy (Vs, solid purple) perturbations. Right panel: Power spectra evolution 
for the pressure (Vs p , dashed red) and nonadiabatic pressure (Vs Puad , solid green) 
perturbations. In both plots, the dashed black line indicates horizon crossing 
corresponding to /cwmap = 0.002 Mpc -1 . 


power spectrum of the curvature perturbation. The amplitude of the entropy 
perturbation power spectrum Vs reduces dramatically when the exchange in 
field dominance occurs at 40 e-folds. It starts to gradually increase until a peak 
is reached at N = 70, at which a drop is then experienced. 


At the end of inflation, we see the magnitude of the entropy power spectrum is 
many orders smaller than the curvature power spectrum. To be specific when 
inflation has ended, the curvature power spectrum is Vq ~ xl0~ 9 whereas the 
entropy power spectrum is Vs ~ 10 -31 . Other final power spectra amplitudes 
are those for the pressure and nonadiabatic pressure perturbations; they are 



r^j 










2.2. RESULTS 


65 


We see that the behaviour of the power spectra evolution of the nonadiabatic 
pressure and entropy perturbation are generally following the same trend from 
horizon crossing, which occurs at N = 15 and onwards. In addition to this, 
we find that Figure 2.2| is in agreement with Figures. 1 and 2 in the works of 
Huston and Christopherson [143J. 


2.2.2 Kinetic coupling 

We have two possible scenarios that will arise in the background dynamics; one 
in which the field 0 reaches the minimum of its potential before the field x, and 
vice versa. We will begin by considering the latter. 


For this scenario, the parameters in the potential are chosen as follows: the 
ratio of the field masses is T = 0.3 with rn^ = 6.395 x 10 -6 . The initial condi¬ 
tions of the two fields are Xini = 12.0 and 0j n i = 11.0 with the field derivatives 
set to slow-roll. In addition, the kinetic coupling between the two fields is 
0 = 0.1. The background dynamics for this set of conditions are displayed in 
Figure 2.3 As shown in Figure [273 , the first 50 e-folds are dominated by the 
field 0 until at which point, there is an exchange in the field contributions, 
leaving the remaining y-field until the end of inflation. 


The evolution of the various power spectra is displayed in Figure |2~4 the left- 
hand panel shows a comparison of the curvature and entropy power spectra 
evolution, and the right-hand panel shows the comparison between the power 
spectra evolution of the pressure and nonadiabatic pressure perturbations. 
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Figure 2.3: Field evolution for the model with one canonical scalar field 0 
(solid red) and one field containing a kinetic coupling x (blue dashed), with 
0 dominating approximately the first 50 e-folds of inflation. The initial field 
values are 0 ini = 11.0 and y ini = 12.0 and masses of the two fields are V = 0.3 
and rritf, = 6.395 x 10 -6 . The kinetic coupling is [3 — 0.1. 




Figure 2.4: Comparison of the power spectra evolution against the number 
of e-foldings for the model with one canonical scalar field and one field with a 
kinetic coupling. For this model, the field 0 dominates the inflationary epoch for 
approximately 50 e-folds. The initial field values are 0i n ; = 11.0 and Xmi = 12.0 
and masses of the two fields are T = 0.3 and m$ = 6.395 x 10“ 6 . The kinetic 
coupling is 0 = 0.1. 

Left panel: Power spectra evolution for the curvature (V-r., blue dashed) and 
entropy (Vs, solid purple) perturbations. Right panel: Power spectra evolution 
for the pressure (Vs p , dashed red) and nonadiabatic pressure (Vs Pnad , solid green) 
perturbations. In both plots, the dashed black line indicates horizon crossing 
corresponding to /cwmap = 0.002 Mpc -1 . 
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The behaviour of the entropy and nonadiabatic component of the pressure per¬ 
turbations is significantly different to the model previously considered in Sec¬ 


tion 2.2.1 However, there is a similarity between the two Figures (2.2 and 2.4). 


For example, a rise and fall is experienced in the evolution of the entropy power 
spectrum, Vs, during the last 6 e-folds before the end of inflation. This feature 
can also be seen in Figure |272j Like the two scalar field model previously stud¬ 
ied, the amplitude of entropy power spectrum at the end of inflation is many 
orders smaller than of the curvature power spectrum. At the end of inflation, 
the final amplitudes for the entropy and nonadiabatic pressure perturbations 
are Vs 10 26 and Vs Pnad ~ 10 50 . Furthermore, the final amplitude of the 
pressure perturbation power spectrum is Vg p ~ 10 -33 . 


We find that the entropy power spectrum amplitude at the end of inflation 
for this model is 10 5 times larger than found in the model considering two 
canonical scalar fields. 


We will now consider the other possible scenario which can arise in this model 
where the field x initially dominates the inflationary epoch. The kinetic cou¬ 
pling acting between the two scalar fields is kept at the same value /3 = 0.1, as 
used in the previous possible set up. However, the parameters used in the poten¬ 
tial are different with the field mass ratio being T = 6.0 and = 1.005 x 10 -6 . 
In addition to this, the starting field values for the two scalar fields are y ini = 7.4 
and 0i n i = 7.5. The background dynamics that are produced by this choice of 
conditions is displayed in Figure 2.5 From Figure 2.5[ it can be seen that the 
final 8 e-folds are dominated by the canonical scalar field. 
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Figure 2.5: Background dynamics for the model with one canonical scalar field 
0 (solid red) and one field containing a kinetic coupling x (blue dashed), with 
X dominating inflation during the first 67 e-foldings. The initial field values 
are 0 ini = 7.5 and Xini = 7.4 and masses of the two fields are V = 6.0 and 
rriff, = 1.005 x 10 -b . The kinetic coupling is 8 — 0.1. 




Figure 2.6: Comparison of the power spectra evolution against e-fold number 
for the case with one canonical scalar field and one field with a kinetic coupling. 
Here, x dominates the inflationary epoch for approximately 67 e-folds. Starting 
field values are 0i n ; = 7.5 and Xini = 7.4 and masses of the two fields are T = 6.0 
and m$ = 1.005 x 10~ 6 whilst the kinetic coupling is (3 = 0.1. 

Left panel: Power spectra evolution of the curvature (TV, blue dashed) and 
entropy (77s, solid purple) perturbations. Right panel: Power spectra evolution 
for the pressure (fPs p , dashed red) and nonadiabatic pressure (Vs Pnad , solid green) 
perturbations. In both plots, the dashed black line indicates horizon crossing 
corresponding to &wmap = 0.002 Mpc -1 . 
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As a consequence of the scalar field containing the kinetic coupling reaching its 
potential minimum earlier during inflation, this has resulted in a one canonical 
scalar field scenario at the end of inflation, which is similar to the two scalar 


fields case studied in Section |2.2.1| Due to this, the shape of the power spec¬ 
tra evolution for all considered quantities will be similar. The power spectra 
evolution for the curvature, entropy, pressure and nonadiabatic pressure per¬ 
turbations are displayed and compared in Figure [2~6j Adding to the statement 
about the similarities in the power spectra evolution between this model and 
that of the two canonical fields, there is a slight difference in the last few e-folds 
in Vs (and Vs Pn!ld ) which relates to the behaviour of the remaining canonical 
scalar field. 


The evolution of the entropy power spectrum is seen initially rising and then 


falling during the remaining e-folds of inflation in Figure 2.2 and Figure 2.4 


Instead we see in Figure |2.6| that the entropy power spectrum continues to 
decrease over time until slow-roll is no longer satisfied. We find the final ampli¬ 
tudes for the entropy and nonadiabatic pressure perturbations are significantly 
larger than for the two standard scalar fields case. However, the values for 
the entropy and nonadiabatic pressure perturbation power spectra amplitudes 
are the similar, if not the same, as those found in the previous scenario, where 
the field x reaches its potential minimum before the field <fi. To be specific, 
the amplitudes for the entropy and nonadiabatic pressure perturbations power 
spectra are Vs r10 27 and Vs PnaA ~ 10 50 . Adding to this, the power spectrum 
amplitude of the pressure perturbation at the end of inflation is V$ p ~ 10 -33 . 
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Furthermore, there appears to be no difference in the curvature and entropy 


power spectra amplitudes between the two scenarios studied here in Section 2.2.2 


2.2.3 A scalar field and DBI field 

We will now study the model containing a canonical scalar field and DBI field. 
In a similar fashion to the previous model containing the kinetic coupling, this 
DBI model will also be associated with the same two scenarios: one in which 
the DBI field decays before the scalar field and the reverse. For all the cases 
considered, the parameters within the DBI model will hold the following values: 
A = 2.0 x 10 12 and // = 0.2 [I53| . 


First we consider the case where inflation is initially dominated by the scalar 
field 0. In this scenario, the field mass ratio is T = 2.0 with the mass of the 
field 0 as = 1.156 x 10~ 5 . The initial values for the two fields are Xim = 2.0 
and 0 ini = 12.0. We display the background dynamics and the boost factor 7 


against the number of e-foldings in Figure 2.7 


In Figure [277) we see that the field 0 reaches its potential minimum when IV = 44 
leaving the DBI field y driving inflation for the remaining e-foldings. From this 
point onwards, the boost factor 7 starts to increase reaching a peak value of 
7 = 540 when N = 59, before it drops slightly to 7 = 490 at the end of inflation. 


In addition to this, the exchange in the field contributions at N = 44 affects the 
power spectra evolution; this is displayed in Figure [ 2 ~ 8 ] through the fall in the 
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Figure 2.7: Background dynamics for the model with a canonical scalar field (ft 
and DBI field y where the scalar field (ft dominates the inflationary period. The 
initial field values are (ft ini = 12.0 and y ini = 2.0 and masses of the two fields are 
T = 2.0 and = 1.156 x 10“ 5 . The parameters required in the DBI model 
are A = 2.0 x 10 12 and fi = 0 . 2 . 

Left panel. Field evolution of the canonical scalar field (ft (solid red) and DBI 
field x (blue dashed). Right panel : Evolution of the boost factor 7 during 
inflation. 


power spectrum amplitude of the entropy perturbation. We see in this model 
that the evolution of the entropy and nonadiabatic pressure power spectra are 
significantly different to the other models previously considered. 


Furthermore, the final value of the entropy power spectrum amplitude, Vs ~ 
10 -16 is markedly greater than found in the two canonical fields model studied 


in Section 2.2.1 Similarly, this increase in amplitude is also found in the power 
spectrum evolution of the nonadiabatic pressure perturbation. For example, in 
this model, the nonadiabatic pressure perturbation power spectrum amplitude 


at the end of inflation is Vs 


Pnad 


10 40 whereas for the two scalar field case 


7V ad ~ lb 54 , which is a difference of 14 orders of magnitude. The power 
spectrum amplitude of the pressure perturbation in this model is Vs p ~ 10 ~ 33 . 
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Figure 2.8: Comparison of the power spectra evolution against the number 
of e-foldings for the model containing a canonical scalar field 0 and a DBf 
field x, where 0 dominates the inflationary epoch. The starting field values 
are 0 ini = 12.0 and Xmi = 2.0 and masses of the two fields are T = 2.0 and 
m ^ = 1.156 x 10 -5 . Parameters needed in the DBf model are A = 2.0 x 10 12 
and /i = 0.2. 

Left panel. Power spectra evolution of the curvature (Vn, blue dashed) and 
entropy (Vs, solid purple) perturbations. Right panel. Power spectra evolution 
of the pressure (Vs p , red dashed) and nonadiabatic pressure (Vs Pnad , solid green) 
perturbations. In both plots, the dashed black line indicates horizon crossing 
corresponding to /cwmap = 0.002 Mpc -1 . 
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Figure 2.9: Background dynamics of the model containing a canonical scalar 
field (j) and DBI field y. In this scenario, the mass of the DBI field x is much 
heavier than the canonical field f , and therefore decays earlier. The initial field 
values are the same as those in the previous case: = 12.0 and y ini = 2 . 0 . 

Masses of the two fields are now T = 35.1 and = 6.50 x 10“ 7 . The parameters 
required in the DBI model are A = 2.0 x 10 12 and p = 0.2. 

Left panel. Field evolution of the canonical scalar field f (solid red) and DBI 
field x (blue dashed). Right panel : Evolution of the boost factor 7 during 
inflation. 


We will now study the other case where the DBI field decays before the scalar 
field. In the potential, the parameters are chosen as T = 35.1 with the mass of 
the 0-ficld as 6.50 x 10“'. The fields will have the same starting values as in the 
previously considered DBI field case. The evolution of the canonical and DBI 


fields alongside the boost factor are presented in Figure |2.9[ and the resulting 
power spectra evolution for the curvature, entropy, pressure and nonadiabatic 


pressure perturbations are displayed in Figure 2.10 


At first, the DBI field dominates the inflationary period until it reaches the 
minimum of the potential well at N = 47, at which point the field will be¬ 
come dominant. The DBI field will now proceed to oscillate about its minimum. 
During these 45 e-foldings of inflation, the boost factor is affected. Specifically, 
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Figure 2.10: Comparison of the power spectra evolution versus e-fold number 
for the model containing a canonical scalar field 0 and a DB1 field x- The DB1 
field x decays earlier during inflation as it is much heavier than the canonical 
field 0. The initial field values are 0 ini = 12.0 and Xini = 2.0 and masses of the 
two fields are T = 35.1 and = 6.50 x 10 - '. Parameters needed in the DB1 
model are A = 2.0 x 10 12 and p = 0.2. 

Left panel. Power spectra evolution of the curvature (Vn, blue dashed) and 
entropy (Vs, solid purple) perturbations. Right panel. Power spectra evolution 
of the pressure (Vs p , red dashed) and nonadiabatic pressure (Vs Pnad , solid green) 
perturbations. In both plots, the dashed black line indicates horizon crossing 
corresponding to /cwmap = 0.002 Mpc -1 . 
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the boost factor increases to a maximum value of 7 = 200 at N = 35 before 


decreasing rapidly to 7 = 1 at N = 47. This is reflected in Figure 2.10 through 
the drop in the power spectrum amplitude of the entropy perturbation 12 e- 
foldings before the end of inflation. 


The behaviour of all power spectra is similar to those seen in the two canonical 


scalar fields in Section 2.2.1 and kinetic coupling models, where the 7 -ficld is 


the first to reach the potential’s minimum in Section 2.2.2 This is due to the 
presence of the remaining field being a canonical field experiencing slow-roll. 
From this, the drop in the entropy power spectrum at N = 47 is due to the DBI 
field reaching the minimum of its potential. The final power spectra amplitudes 
for the entropy and nonadiabatic pressure perturbations are Vs rs./ 10 37 and 
V$ p nad ~ 10~ 62 . For completeness, the pressure perturbation power spectrum 
amplitude is Vs p rv_/ 10 " 34 . 


The final amplitudes for the entropy and nonadiabatic pressure perturbations, 
Vs and Vg Pjmd for this scenario are much smaller than in the two scalar field 
case in Section 12.2.11 


We find in the model containing the kinetic coupling that regardless of which 
field reaches the minimum of its potential first, it has no effect on the power 
spectra generated. However, in the DBI model considered in this section, we 
find that the ordering of field decay has an effect on the power spectra produced, 
in particular, the power spectra of the nonadiabatic pressure perturbation and 
ultimately, the entropy perturbation. We find that the difference between nona- 
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diabatic pressure power spectrum amplitude of the two scenarios considered here 
is incredibly large. To be precise, in the scenario where the canonical field is 
the first to reach the potential minimum, the power spectrum amplitude of the 
nonadiabatic pressure perturbation is Vs Pnad ~ 10 -40 , and for the alternative 
scenario, the amplitude is Vs Pnild ~ 1CT 62 . Naturally, the amplitude of the en¬ 
tropy power spectrum is also greatly affected; when the canonical field reaches 
the minimum first, the entropy power spectrum amplitude is Vs rv_/ 10 16 , and 
vice versa, the amplitude is Vs ~ 10 -37 . There is no difference in the pressure 
perturbation power spectrum amplitude between the two DBI scenarios. 


2.3 Conclusion 

We have studied the evolution of the nonadiabatic pressure perturbation <5p nad in 
the context of multi-field inflation, furthering the work of Huston and Christo- 
pherson [1143} . In this work, we have examined three models where each model is 
composed of two massive scalar fields. However, they all differ in kinetic terms: 
first model consisting on two canonical scalar fields, second model with one 
canonical scalar field and one field with the presence of a kinetic coupling, and 
finally, the third model containing a canonical scalar field and a DBI field. For 
all models considered, inflation is driven by both fields until one field reaches 
the minimum of its potential, at which point, inflation is then driven solely by 
the remaining field. 

In the process of this numerical study, we confirm our results for the model 
composed of two canonical scalar fields are in agreement with Huston and 
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Christopherson jl43j . In addition to this, we found that the nature of the 
kinetic term affects the evolution and ultimately, the final value of the nona- 
diabatic pressure perturbation. For example, in the case involving the kinetic 
coupling, we find that the final amplitude of Vs Pnad is independent of which field 
primarily decays. Yet, an increase in the kinetic coupling generates an increase 
in V$ PnaA . In particular, the difference between no and a nonzero, specifically 
when f3 = 0.1, kinetic coupling yields an increase in the final amplitude of Vs Pnad 
of roughly five orders of magnitude. 

The background dynamics, specifically the speed at which the first field reaches 
the minimum of its potential, greatly influences the evolution of Vs- At the 
point where the first field reaches its potential minimum, it corresponds to an 
increase in Vs, and in turn affects V-r. We show that this effect is most ap¬ 
parent when considering the case of two canonical scalar fields. For all other 
models considered here, the decay of the first field is not as dramatic as that of 
the canonical case, and hence, the changes in Vs are less visible. 

In the case of the canonical scalar and DBI fields, we see that if the DBI field 
drives the remaining e-folds of inflation, it causes the final amplitude of Vs Pnad 
to be significantly larger than other models examined. On the contrary, if the 
DBI field has reached its potential minimum earlier on during the inflationary 
epoch, and therefore is not significant for last few e-folds of inflation, we see that 
Vs p nad is similar to the canonical case. Furthermore, we show that the DBI field 
driving the last few inflationary e-foldings produces the largest amplitude of Vs- 
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CHAPTER 2. NONCANONICAL MULTIFIELD INFLATION 


The nonadiabatic pressure perturbation can be shown to source vorticity at 
second order in perturbation theory, and in turn will impact CMB B-mode 
polarization predictions. We have shown that a model containing a DBI field 
driving the remaining e-folds of inflation will provide the largest final amplitude 
for ’P« 5 Pnad , and hence supply a larger vorticity source when compared to other 
models considered here. Further study into the mechanism of both pre- and re¬ 
heating in these models is imperative in order to estimate the amount of entropy 
perturbations in the radiation dominated epoch; from [156] . it was shown that 
in this epoch, the nonadiabatic pressure perturbation can be generated. The re¬ 
sults will depend on the model considered due to the particular features within, 
for instance on couplings between the inflaton ficld(s) and matter fields. 



Chapter 3 

Implications of sharp transitions 
on the effective Planck mass 
during inflation 


Cosmological inflation has been widely accepted as the process which gives rise 
to the formation of structure in the Universe. However, there are many in¬ 
flationary models which satisfy the vast quantity of observational data that 
have been obtained through experiments. This makes determining the defini¬ 
tive model of inflation incredibly difficult. One area of study which has been 
active are scalar-tensor theories. 

Scalar-tensor theories are the most natural extensions to Einstein’s Theory of 
General Relativity and were first studied by Jordan in 1959 |157j . In 1961, 
Brans and Dicke published a paper suggesting an alternative to Einstein’s The¬ 
ory of General Relativity [102] as a way to describe classical gravity. Within 
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this theory, there are gravity-matter couplings and from this, two frames arise: 
the Jordan frame and the Einstein frame. The Jordan frame is a frame in which 
the Einstein-Hilbert action has been modified to include a coupling to a scalar 
held. By performing a conformal transformation of the metric [ 158 . 159] . the 
Jordan frame can be related to the Einstein frame. In the Einstein frame, the 
Einstein-Hilbert action is as usual with the additional of a scalar held. The 


action written in these two frames can be seen in Chapter 1 in Section 1.6 enti¬ 


tled “Alternative models”: Eq. (1.99) for the Jordan frame and Eq. (1.104) for 
the Einstein frame. Scalar-tensor theories provide interesting phenomenological 
models; they can be used to study the variation of Newton’s constant G. 


In 1991, Accetta and Steinhardt stated that if the gravitational constant os¬ 
cillates at high frequencies when compared to the Hubble expansion rate, it can 
affect cosmological measurements [160]. The authors state the only requirement 
is that the system contains a massive scalar held which is coupled nonminimally 
to gravity. Oscillations in the gravitational coupling are obtained through the 
scalar held oscillating about its ground state value. The frequency of the scalar 
held oscillations are related to the mass of the held. In their 1991 paper, the 
authors focus on epochs after inflation, i.e. the radiation and matter dominated 
epochs. This work was extended in 1994 by Steinhardt and Will who considered 
the variation of the the gravitational coupling after inflation [ 161 ]. The authors 
constructed a model with oscillations in Newton’s constant by using a Brans- 
Dicke model with a nonminimally coupled massive scalar held 0, which has 
the ability to be displaced from the minimum of its effective potential during 
inflation. Inflation in this model is not driven by the held 0, but instead by an 
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inflaton field within the matter sector. During the inflationary epoch, the cou¬ 
pling between the gravity and scalar field sectors generates a finite displacement 
of the held 0 from its potential minimum, creating a potential minimum which 
has a false vacuum energy density. Once inflation has ended, reheating starts, 
which causes the inflaton energy density to be converted into radiation. Due to 
this, the potential minimum now occurs at 0 = 0 resulting in the held having 
two options: to oscillate or to slowly roll. The choice is determined whether 
the period of oscillation is greater than the expansion time scale. For the case 
when the period of oscillation is greater than the expansion time scale, the held 
begins by slow-rolling down the potential eventually leading to oscillations. For 
the reverse case, the held immediately starts with oscillatory behaviour. 

It is known that the gravitational constant could not have varied too much 
during the Universe’s lifetime. In fact, the gravitational constant could have 
only varied by [ 1621 l63 j 

[U < 10- 12 yr-‘ . (3.1) 

Or at 

There are many experiments that measure the time variation of the gravita¬ 
tional constant: lunar ranging observations [iMi. rm. Big Bang Nucleosynthesis 
(BBN) [ 166Ul67I'168l 169], gravitational waves |170j . and more recently WiggleZ 
m ■ Other experiments that constraint the varying Newton’s constant involve 
stellar objects; some examples include comparing the ages of individual Globu¬ 
lar Clusters to the age of the Universe m, observations of Type la supernovae 
|17311174] . pulsating white dwarfs |17511176] . pulsars |17711178111791118U] and 
neutron star surface temperatures [181], 
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As we will see in this chapter, varying the gravitational coupling can affect the 
scalar-to-tensor ratio. In the first quarter of 2014, the BICEP2 team reported 
a detection of B-mode polarization of the CMB; this was published later in 
2014 cm When this detection is interpreted as gravitational waves produced 
during the inflationary epoch, it yields a tensor-to-scalar ratio r of 

r = 0.2(tf“;»; , (3.2) 

at around I ~ 80 which corresponds to a scale of &bicep 2 — 0.005 Mpc^ 1 . 
However, this is in contention with the Planck results from 2013 ra stating 
an upper limit of r < 0.11 (95% C.L.) at the scale fcpianck — 0.002 Mpc^ 1 i.e. 
I ~ 28. These values suggest that the scale of inflation is at the GUT scale and 
furthermore, a tensor-to-scalar ratio that is r > 0.11 will rule out a significant 
number of inflationary models such as those where the field displacement is less 
than the Planck mass A0 <g; Mpi P21HS31 : models where r > 0.1 is obtainable 
include those inspired by string theory: assisted inflation [TMi. 185 . 11861 [187] . 
N-flation [188, 189] . a model similar to the latter [1901 191, 192] and monodromy 

(Mi iisa iffl. 

Since the announcement of the BICEP2 results, many authors have proposed 
solutions to resolve the discrepancies between these results and from Planck 
on both theoretical |19611197] 1198] 1199] 1200112011120211203] and experimental 
grounds [204] [US |206j. With the upcoming Planck 2014 results and those 
from experiments scheduled for years to come, they will continue to test exist¬ 
ing models of inflation; other searches for features in the power spectra include 
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[8811891 1207112081120911210] . 


In this chapter, we consider an action of scalar-tensor form, specifically study¬ 
ing the effects of a sudden transition in Newton’s gravitational constant during 
inflation and its consequences on the scalar curvature and tensor perturbations. 
The transition made is smooth step in nature and not violent, and so the vari¬ 
ations of Mpi are not of order one, but rather of order a percent or less. Two 
scenarios will be considered in this chapter: one in which the Brans-Dicke field 
will play the role of the inflaton, and a two-field model where this role will 
be assigned to the second auxiliary field. We begin by presenting the action, 


background and perturbation equations for the single-field case in Section 3.1 


and the model in Section 3.2 which contains three benchmark models in Sec¬ 


tion 3.2.1 The addition of a minimally-coupled field to the action is considered 


in Section 3.2.2 In this section, the relevant equations are presented alongside 


one case. We conclude with Section 3.3 containing a summary of our findings. 


3.1 Field equations 


3.1.1 The background 

The action that we consider is of the form 



d *Xyf=Q 


Ml i 


h<p)r 




(3.3) 


where Mp\ is the bare reduced Planck mass with Mpi 2 = (87 tG*) 1 , and con¬ 
tained is the bare gravitational constant G*, R is the Ricci scalar and U(ip) is 
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the potential of the scalar held p . The gravitational sector is conpled to the 
scalar held via the effective Planck mass F(tp), i.e. Mp j x F{tp). From now on, 
the bare reduced Planck mass is set to one Mp\ = 1. 


All calculations will be performed in this frame known as the Jordan frame; 
this will allow us to model the variations in the Planck mass with ease. Physi¬ 
cal quantities such as the power spectra are the same in both the Jordan frame 
and the Einstein frame (see [2TT[ : 212 , 213] and the references therein). 


By varying the action in Eq. (3.3) with respect to the metric g^ v yields the 
Einstein equations 


F(ip) L 


d v p - -g liv g a0 d a p d 8 p 
+ V^V u F(<p) - g^UF{p>) - g^ v U(p>) 


(3.4) 


where G IJiy is the Einstein tensor and □ = g^V^Nu is the covariant d’Alembertian 
operator. 


The equation of motion of a scalar held is obtained by varying the action in 


Eq. (3.3) with respect to such held; for the held ip : it is as follows 


Of = Ujf(tp) - ^F iV (f)R , 


(3.5) 


with F = dF/dtp, noting that denotes partial derivatives with respect to the 
held (p. This background equation can be rewritten in a different form known 
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as the Brans-Dicke equation. We first require the Ricci scalar 11. By taking the 


trace of the Einstein held equations in Eq. (3.4) and using the definition of the 


Einstein tensor as stated in Eq. (1.5), the Ricci scalar is 


R = 


F(V) 


d^ip + 3 +4U (if) 


(3.6) 


with 


DF(ifi) = F iV (if)nip + F^^g^a^ipd^ip . 


(3.7) 


Substituting this for R in Eq. (3.5) yields the Brans-Dicke equation 


2wUp = -w^g^d^ipd v p - 4F iip (ip)U (ip) + 2 F(ip)U tip (ip) , (3.8) 


where w = F(<p) + | F yip {(pf 


We assume that the Universe is both homogeneous and isotropic on large scales, 


and therefore we use the FRW line element as stated in Eq. (1.8) in Section 1.3 


Using this, the Brans-Dicke equation in Eq. (3.8) is 


ip + 3 Hip = — l-w )ip ip 2 + 4 F^U - 2 FU V \ . 
2w 


(3.9) 


where dots represent derivatives with respect to cosmic time. Furthermore, 
due to the form of the Einstein held equations, the standard Friedmann and 


acceleration equations are obeyed — see Eqs. (1.9) in Section 1.3 — with the 
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effective energy density and effective pressure given by 



l -V 2 + U-2>HF , 
- U + F + 2 HF 


(3.10a) 

(3.10b) 


It is important to note that these are effective quantities and that the corre¬ 
sponding energy-momentum tensor Tj£/ behaves as a perfect fluid and hence 
obeys = 0, and the leading diagonal is T ,J V = diag(— p,p,p,p), where p 

and p denotes the energy density and pressure of the perfect fluid. 


In order to extend and test the generalities of this single-field theory, we later 
consider a two-field model. The two-field model will contain a second minimally- 
coupled field x in conjunction to the single-field case. The action of the auxiliary 
field is 


g(x) 



g^d^x duX + 2V(x) , 


(3.11) 


with V(x) denoting its potential, which is of the quadratic form: 


v ( x) = § m h 2 


(3.12) 


We shall first proceed with the single-field model and later in the chapter, 


specifically in Section 3.2.2 introduce the second scalar field into the model. 
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3.1.2 Perturbations 


We will now focus our attention to the first-order perturbation equations, which 


will be studied in the Newtonian gauge as presented in Section 1.5.2 For 
convenience, we now remind the reader of the line element in the Newtonian 
gauge as first stated in Chapter 1 


ds 2 = —(1 + 2'h)d t 2 E a(f) 2 (l — 2E)S i jdx 1 'dx j , 


(3.13) 


where and <f> are the scalar metric perturbations. In addition, we will use the 
held decomposition, where we split the held into its background and perturbed 
parts, and work with Fourier modes. 


By perturbing the metric, the perturbation equations of the two fields are 


Sip E 


3H E ’—Al 
w 


E 


Sip 


w 


1/1 


E 

,<p 


r +TS 


VO 


k 2 


— (\k E 3 < h)0 E — 

ZZ 7 


- 4 FJJ E 2 FU 




E = 0 , 


Sip 


(3.14) 


In the Newtonian gauge, the perturbed Einstein equations are given by the 
following 


k 2 1 

3 H{EEHE)E^E = --Sp , 
a z 2 

(3.15a) 

EE HE = ~Sq , 

(3.15b) 

EE(2HE 3 H 2 )E E H(E E 3$) = ^Sp . 

(3.15c) 











88 CHAPTER 3. TRANSITIONS IN THE EFFECTIVE PLANCK MASS 


The right-hand side of these equations are related to perturbations of the 
effective energy-momentum tensor T^, and are equal to perturbations of the 
energy density, momentum potential and pressure, respectively 


5p 


5q 

5p 


1 

F 


<p5(p - 0 2 T + UJ(p + 3F(<1> + 2HH) - 3 H{SF + H5F ) 


k 2 

- 2 5F 
a 2 


1 

F 


pSp + 5F — FT — H5F 


(3.16a) 

(3.16b) 


c pS(p - ip 2 ^! - U v 5(p — pSF + 5F + 2 H5F - M - 2Fd> 


u 2 

-2(F + 2TTF)T + —SF 
a 2 


(3.16c) 


By using the ij-component of the Einstein equations, we find that anisotropic 
stress is present in the Jordan frame 


F F 


(3.17) 


Observables are required to relate theories to experiments and examples of such 
observables include the spectral index n s and its running a. For convenience we 
will now recall parts of Section [L5[2 entitled “Scalar modes”. Both the spectral 
and running indices are related to the scalar perturbations through the power 


law, which we recall from Eq. (1.76) 




\ n s (ko) — 1+1 ln(fe/feo) a 


(3.18) 
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where k 0 is a pivot scale and V^(k 0 ) is the amplitude of the scalar power spec¬ 
trum at such pivot point. 


The scalar perturbations are provided by the curvature perturbation on con¬ 


stant hypersurfaces ( as defined in Eq. (1.59) 


C = - —V , 

P 


(3.19) 


and the resulting power spectra generated is 


P 


(3.20) 


The pivot scale chosen in this work is that used by the Planck experiment which 
is k 0 = 0.002 Mpc -1 , and the amplitude of the scalar curvature power spectrum 
at this pivot scale is *P^(/cq) ^ 2.15 x 10 '* m- Furthermore, the running index 


a is defined as in Eq. (1.77) 


dn s 


a = 


(3.21) 


d In k 

From the Planck experiment [214] . the current best fit values for both the 
spectral index and its running are 


n s = 0.9603 ± 0.0073 , a = - 0.0134 ± 0.0090 


(3.22) 


In addition to considering scalar perturbations, we shall also study tensor per¬ 


turbations. 
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We will briefly recapitulate material on tensor perturbations from Section 1.5.4 
for practical purposes. The equation of motion for tensor perturbations is given 


in Eq. (1.92) 


h" + 2-ti + k 2 h = 0 , 
a 


(3.23) 


with the tensor power spectra given by 


U3 

Pt = 2^| 2 


(3.24) 


Both these equations will have the standard oscillatory initial conditions 


The tensor-to-scalar ratio r is used to relate the tensor and scalar power spectra. 
It is defined as follows da 


r = 


8 Vt 


(3.25) 


3.2 Model with step variation in the Planck 
mass 

We will now examine the smooth step transition in the effective Planck mass 
_F(<p); the effective Planck mass and the potential considered are 


F((p) — 1.0 — /3 { 1.0 + tanh[( cp — V 7 *)/^]} > (3.26) 

U(<p) — ^rn^ip 2 , (3.27) 


where is the mass of the scalar field, (3 is a dimensionless constant, and both 
7 and <p* are constants of mass dimension. The parameters /3 and 7 determines 
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the sharpness and amplitude of the transition and the parameter 9 ?* determines 
the field value at which the transition occurs. We have chosen the quadratic 
potential for concreteness, but we expect the features observed due to the sharp 
transition to also appear in other choices of potential. 

In the next section, we will firstly consider a single-field model, in which the 
Brans-Dicke field drives inflation. Following this, we will then consider a two- 
field model where we introduce a second minimally-coupled scalar field which 
will act as the inflaton. For the two cases and their benchmark models studied 
here, the values of the parameters are chosen so that successful inflation is ob¬ 
tained with the inflationary period lasting 66 e-folds. 


Before proceeding, we will now show that the Jordan frame model that we 
are studying is not equivalent to the Einstein frame model of an inflaton po¬ 
tential with step as described in ||87jj. Through the conformal transformation 


as outlined in Section 1.6, we can use this to transform the model considered 


in Eqs. (3.26) and (3.27) to the Einstein frame; the potential in the Einstein 


frame is 


u ( 0 ) = u m = _ m2 y _ 

F 2 (<p) 2 {1.0 — /3 { 1.0 + tanh[(</p - <~P*)h]}} 2 


(3.28) 


It is important to note that the potential above and that used in [57] are the 
same when /3 -C 1. However, there will be different dynamics due to the 
differences in the kinetic terms. When the model considered in this chapter 
is transformed to the Einstein frame, it will have noncanonical kinetic terms 
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unlike that studied in [87j. The canonical field p is related to the noncanonical 


field p via the relation — see Eq. (1.101) 


'2F{py + 3F^(py 

1 »'- 2 F(„f iV 


(3.29) 


By rewriting the equation above, the field p can be expressed in terms of p 
at least implicitly. The potential U(p) can now be rewritten in terms of the 
canonical field p. However, the form of the potential will not longer be a step 


potential as that stated in Eq. (3.28). As a result, the Einstein frame model 


containing a step potential is not the equivalent to a Jordan frame model con¬ 
taining a step transition in the effective Planck mass; the model considered in 
this chapter is different to those previously studied in the literature, and hence 
the results from this study will provide significantly different predictions for the 
scalar and tensor power spectra. 


The dynamics of the fields will be numerically solved, following the method 
outlined in EH; the derivative of the background fields are given their slow- 
roll values, and the initial field perturbations will have the standard oscillatory 
Bunch-Davies initial conditions [87]. In order to calculate the tensor perturba¬ 
tions generated by the system, we employ the methods described in [SZl 'HI. 

For the effects of this type of transition to be observable, they must occur be¬ 
tween 60 and 50 e-folds before the end of inflation [ 154. ' 218j . This will translate 
to a range of scales observable using experiments such as Planck. These step 
changes in the Planck mass could result from a first-order phase transition in 
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the vacuum expectation value (VEV) of a Brans-Dicke field [219 ]. A similar step 
transition in the inflaton’s potential was considered in [87] (see [220 ] for a sharp 
transition in the effective Planck mass) and was shown to result in oscillation 
in the primordial scalar curvature power spectrum. However, no oscillations in 
the scalar power spectrum were seen when considering a nonminimally coupled 
Brans-Dicke field with a transition in the potential. 


3.2.1 Single-field model 

Minimally-coupled limit 


We shall begin with the minimally-coupled case; in this case, the parameter /3 
is zero and hence F — 1. The starting value for the ip field is <p ini = 16.179 with 
mass m,_p = 6.5 x 10 -6 . The mass is chosen so that the power spectrum for the 
scalar perturbations at the Planck pivot scale is approximately 2.15 x 10~ 9 . The 


power spectra for the scalar and tensor perturbations are given in Figure 3.1 


Notice that we have defined the number of e-folds N such that N = 0 at the 
start of inflation. With the coupling between the gravitational sector and scalar 
field set to one, there are no features generated in this model, as expected. For 
this case, the spectral and running indices are calculated to be 


n s = 0.968865 , a = 0.00107427 . (3.30) 


In addition to this, the tensor-to-scalar ratio at the Planck pivot scale is 


r(fc 0 ) = 0.133205 . 


(3.31) 
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Figure 3.1: The power spectra ( left panel ) of the scalar (T% solid black) and 
tensor (Vt, blue dashed) perturbations and the associated tensor-to-scalar ratio 
(right panel) against wavenumber for the minimally-coupled model, i.e. (3 = 0. 
The parameters used are the initial field value, </q n i = 16.179 and field mass 
= 6.5 x 10 -6 . 


For this case, a quadratic potential was chosen for simplicity. The model is 
slightly under pressure from the Planck experiment [214] (cf. 221] as an attempt 
to reconcile the model with the Planck data), although it is still within the 68 % 
C.L. in the n s — r plane. 


Benchmark 1 


In the first of three benchmark models, we will consider a steep transition in the 
effective Planck mass, which is reflected in a violent feature in the slow-roll pa¬ 
rameter e; this evolution in the slow-roll parameter and of the effective Planck 


mass against the number of e-folds is displayed in Figure 3.2 The parameters 


chosen are (3 = 0.0460, 7 = 0.145 and (p* = 15.8, with the mass of the field as 
nip = 2.1 x 10 -5 and the initial field value is p- m \ = 16.5783. 


The resulting scalar and tensor power spectra are displayed in the left panel 
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of Figure |3.3| We see an extremely sharp dip in the power spectrum of the 
scalar perturbation at k ~ 0.003 Mpc -1 , which is followed by damped oscilla¬ 
tory behaviour for smaller wavelengths. For this set of parameters, there is a 
noticeable feature in the tensor power spectrum, as demonstrated in Figure [3l4| 
Fading oscillatory features in the primordial scalar power spectrum can also 
occur from jumps in the potential [EH [2221 [223, 223 , 1225] . particle production 
during inflation [ 2261 227 . 1228 ] or turns in the inflaton trajectory in the land¬ 
scape of heavy fields [222, 230 . [23Ti,1232] . 


The clear differences between the scalar and tensor power spectra in Figure 3.3 


can be explained by examining the behaviour of the slow-roll parameter e. For 
the single-held case only, the scalar power spectrum follows as ~ H 2 /e, 
whereas the tensor power spectrum follows as V t ~ H 2 . As a result, a sharp 
dip in the scalar power spectrum is observed at the same time when a sharp 
peak occurs in the slow-roll parameter, however, the tensor power spectrum is 
largely unaffected. 


Using both power spectra, the tensor-to-scalar ratio can be calculated and its 
behaviour as a function of the wavenumber is presented in the right-hand side 


panel of Figure 3.3 Although a very large unrealistic tensor-to-scalar ratio is 


generated in the region of the Planck pivot scale, this can be useful in constrain¬ 
ing the parameters in the nonminimal coupling F(ip). 
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Figure 3.2: The evolution of the slow-roll parameter e (left panel ) and the 
nonminimal coupling F{<p) (right panel ) for the first 20 e-folds of inflation for 
benchmark 1. The initial field value is <p- m \ = 16.5783 and mass is m v = 
2.1 x 10 -5 . Model parameters required in the nonminimal coupling F(ip) are 
(3 = 0.046, 7 = 0.145 and </?* = 15.8. 


Benchmark 2 


In this model we show that by selecting the appropriate model parameters, it 
is possible to find an agreement between the results from the Planck and BI- 
CEP2 experiments. The Planck experiment has placed an upper bound on the 
tensor-to-scalar ratio of r < 0.11 at the pivot scale k 0 = 0.002 Mpc -1 , whereas 
BICEP2 has stated a value of r = 0.2 at the scale /cbicep 2 = 0.0048 Mpc -1 . The 
two observations would not agree with the standard power law power spectrum. 
However, through the transition in the effective Planck mass, the disagreement 
between the two experiments can be weakened. The parameters for such a 
model are = 6.9 x 10 -6 , f3 = 0.002, 7 = 0.111 and = 15.49, with the 
initial field value = 16.2271; the gravitational coupling F((p) and the evo¬ 
lution of the slow-roll parameter e is displayed in Figure 3.5| 


In Figure 3.5, we see that the slow-roll parameter creates a peak due to the 
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Figure 3.3: The power spectra (left panel ) for the scalar (7^, solid black) and 
tensor (Vt, blue dashed) perturbations and the associated tensor-to-scalar ratio 
r (right panel ) against wavenumber, for benchmark 1. The initial field value 
is (^i n i = 16.5783 and field mass is = 2.1 x 10 -5 . Parameters used in the 
nonminimal coupling F(ip) are [3 = 0.046, 7 = 0.145 and 99 * = 15.8. 


increase in the coupling at approximately N — 6 . As previously seen in bench¬ 
mark 1 , a small effect on the slow-roll parameter e greatly affects the tensor-to- 
scalar ratio; notice that, with this combination of parameters, the initial value 
of F(</?) = 0.996, deviating by less that 0.5% from minimal coupling, compared 
with ~ 10% deviation in benchmark 1. The resulting power spectra for the 


scalar and tensor perturbations are displayed in the left panel of Figure 3T and 
using the two power spectra, the tensor-to-scalar ratio can be calculated; the 
tensor-to-scalar ratio versus wavenumber is shown in the right panel of Fig¬ 
ure 


3.6 As expected, a reduction in the power of the scalar power spectrum 


is observed whereas the tensor power spectrum is unaffected. In Figure 3.6 


we have demonstrated that the tensor-to-scalar ratio can be restricted to the 
allowed values constrained by BICEP2. The maximum value of the spectral 
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k [MpcT 1 ] 

Figure 3.4: The tensor power spectrum Vt of benchmark 1. The starting field 
value is p > ini = 16.5783 and field mass is m v = 2.1 x 10 -5 and in the coupling 
F((p) are (3 = 0.046, 7 = 0.145 and tp* = 15.8. 


index in the vicinity of the feature is obtained by means of [233 


n s (k) 


~ 1 + 


din V c (k) 
d In k 


(3.32) 


and for k > 10~ 2 Mpc _1 , the spectral index is n s ~ 0.98. One can obtain a 
rough estimate on the maximum magnitude for the nonlinearity parameter /nl 
in the squeezed limit from that [ 234 ] 


|/NL|max — -^2^ 


n s 


~ 0.45 . 


(3.33) 


The dimensionless parameter /nl measures the amplitude of non-Gaussianity 
- see the following references: [ 235 , 12361123711238] . Thus, even in the vicinity 
of the transition, this model remains consistent with the Planck limit for local 
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Figure 3.5: Background dynamics for benchmark 2 model: the evolution of 
the slow-roll parameter e ( left panel) and the corresponding F(tp) (right panel) 
for the first 20 e-folds of inflation. The initial held value is <^j n i = 16.5783 with 
a mass of m v = 2.1 x 10 -5 and the parameters required in the nonminimal 
coupling F(<p) are f3 = 0.046, 7 = 0.145 and 99 * = 15.8. 


non-Gaussianity of = 2.7 ± 5.8 |239j . 

Benchmark 3 

Finally, we have also studied the case when the factor /3 within the gravita¬ 
tional coupling F(cp) is of the opposite sign i.e. /3 is negative. For this bench¬ 
mark model, the parameters used in the nonminimal coupling are (3 = — 0.005, 
7 = 0.100 and ip* = 14.64. In addition, the starting held value is = 15.9055 
and mass used is = 6.5 x 10~ 6 . Both the evolution of the slow-roll parame¬ 
ter £ and coupling F(ip) for this choice of parameters are displayed in Figure [T7] 


From Figure |3.7 we see the change in sign in /3 causes a dip in the evolu¬ 
tion of the slow-roll parameter, and from this, we expect a bump in the scalar 
power spectrum; the power spectra for both the scalar and tensor perturba¬ 
tions, and the tensor-to-scalar ratio with respect to wavenumber are presented 
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Figure 3.6: The power spectra (left panel ) for the scalar (V^, solid black) and 
tensor (Vt, blue dashed) perturbations and the associated tensor-to-scalar ratio 
r (right panel ) against wavenumber for benchmark 2. Parameters used in the 
nonminimal coupling are f3 = 0.002, 7 = 0.111 and 92 * = 15.49. The starting 
field value is = 16.2271 and field mass is rn^ = 6.9 x 10 -6 . 


in Figure 3. 


We see the sign change in [3 accompanied by this choice in model 


parameters causes the dip to last for approximately 7 e-folds. 


The reduction in the slow-roll parameter has had no considerable effect on the 
tensor power spectrum. Furthermore, this reduction is the source of enhance¬ 
ment in the scalar power spectrum in the range: 10 ~ 2 Mpc -1 < k < 10 1 Mpc -1 , 
which in turn suppresses the tensor-to-scalar ratio in the wavenumber range: 
10 -1 Mpc -1 < k < 10 1 Mpc -1 . 
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Figure 3.7: The evolution of the slow-roll parameter e (left panel ) and nonmin- 
imal coupling F{<p>) (right panel) for the first 20 e-folds of inflation for bench¬ 
mark 3. In this case, the initial field value is = 15.9055 and mass is 
rritp = 6.5 x 10 -6 . The parameters required in F(<p) are B = — 0.005, 7 = 0.100 
and 99 * = 14.64. 


3.2.2 Two-field model 


To extend the single-field model given by the action in Eq. (3.3), a second 


minimally-coupled field as stated in Eq. (3.11) is added. This addition was 


made in order to test the results generated by the single-field model, i.e. the 
distinctive feature appearing in the primordial scalar power spectrum as seen 


in Figure. |3.6[ for example. In this two-field model, the step transition again 
occurs in the Brans-Dicke field 99 , however the inflaton is played by the auxiliary 
field x- The Einstein field equations takes the form 


- F 


T !$+dn<p d„<p--g^g al3 d a ip ^V^F-g^UF-g^U 


, (3.34) 


where 


t(x) _ 

jlV 


2 5SG) 


V~~9 s g lw 


(3.35) 
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Figure 3.8: The power spectra (left panel ) for the scalar (V^, solid black) and 
tensor Vt, blue dashed) perturbations and the associated tensor-to-scalar ratio r 
against wavenumber, for benchmark 3. Model parameters used are d — — 0.005, 
7 = 0.100 and 99 * = 14.64. The starting field value is v?i n i = 15.9055 and field 
mass is rn v = 6.5 x 10~ 6 . 


is the energy-momentum tensor of the field y. 


Varying the full action (a combination of Eqs. (3.3) and (3.11)) with respect to 
the two scalar fields yields their equations of motion as shown below 


2 wUip = F^T™ - - AF tV U + 2 FU V , (3.36a) 

□X = v ,x • (3.36b) 

Applying the FRW metric to these equations produces the background equa¬ 
tions of motion 


- 2 w{Cp + 3 Hip) = F^(\ 2 - 4V) - zu itp (p 2 ~ ^F jip U + 2FU iV , 
X + 3 Hx + V x = 0 . 


(3.37a) 

(3.37b) 
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The Friedmann equations take the following forms: 


H 2 = — 
3 F 

-V 2 +-X 2 + U+V-3HF 


p> 2 + x 2 + F-HF 



(3.38a) 

(3.38b) 


The relevant perturbation equations may be found in 2 10 . 12T1J and are given 
by 


Sip + 


LU / 


3 H H- —ip Sip 


w 


jdx) + I ( 


w 


(T T Sd*)^ T — 

G7 


W 


1 ( 1 


^--A-^U~2FU v ) 


F m tM - AF_JJ + 2 FU, 




+ = 0 , 

ZTZJ 

Sx + 3 HSx + k J\ - (T + 3<f>)x + 2E X T + V xx S X = 0 , 


k 21 

H —2 
a 2 


Sip 


(3.39a) 

(3.39b) 


where T^ is the trace of the energy-momentum tensor of the field x and 8TW 
is its perturbation. The right-hand side of the perturbed Einstein equations in 
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Figure 3.9: Evolution of ip and x versus e-fold number N (left panel), and the 
field trajectory in the tp — x plane for the same e-foldings, illustrating the sharp 
turn at IV ~ 10 (right panel). 


Eq. (3.15a) are given by 


6 P= E 


_ J-q, + - x 2 * + (U iV 5<p + VJx) + 3F(<h + 2 HU) 


k 2 

3 H{5F + HSF) - — 5F 


Sq = 
5p = 


1 

F 


<p8<p + x8 x + 8F-F^ - H8F 


(3.40a) 

(3.40b) 


<p5<p -Jq, + tfx - x 2 ^ - (U >V 6<P + VJx) - PtotSF + 8F + 2 HSF 


k 2 

F'F - 2- 2 (F + 2 HFJ + —8F 


(3.40c) 


where p to t is the total effective pressure from the two fields. It is defined as 


Ptot 


F 


1 -2 , 1 • 2 

2 V + 2 X 


U-V + F + 2 HF 


(3.41) 


For this model the parameters in the coupling F{<p) are chosen as f3 = 0.009, 
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7 = 0.111 and 92 * = 1.0. The initial field values are = 4.5 and x ini = 15.489, 


and their respectively masses are = 4.48 x 10 5 and m x = 5.6 x 10 b . As 
a result, the background dynamics of the two fields ip and x are displayed in 


Figure [379] and it is clear that the final 50 e-folds of inflation are driven by the 
field x■ The evolution of the slow-roll parameter e and the effective Planck mass 


F are shown in Figure 3.10 


In Figure 3.10 we see the feature from the sharp transition in F(ip) appears 


upon the rolling of the Brans-Dicke field towards the origin. This occurs at 
approximately 7 e-folds from the start of inflation. We see that the fluctuation 
in the slow-roll parameter e is of comparable nature to that displayed in bench¬ 
mark 1 . 


We observe from Figure 3.11 a smooth enhancement in the tensor power spec¬ 
trum for longer wavelengths i.e. k < 10 _ 2 Mpc _1 . This can be understood 
by the overall reduction in the slow-roll parameter after the departure of the 
second field. The ^-dependent tilt of the tensor power spectrum is given by 
[721 12421 12431J 


n t = -2e - 


HF 


(3.42) 


The gradient in the nonminimal coupling F(tp) is zero before and after the 
transition, therefore, the slope of the tensor power spectrum is provided solely by 
the slow-roll parameter e. This parameter is considerably larger pre-transition. 
In contrast to the single-field model, we do not see the effect of the sharp feature 
in the effective Planck mass appearing in the scalar power spectrum for the two- 
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Figure 3.10: Evolution of the slow-roll parameter £ (left) and the effective 
Planck mass F (right) for the first 20 e-folds of inflation in the two-held model. 
The model parameters are m x = 5.6 x 10 -6 , m v = 4.48 x 10 -5 , /3 = 0.009, 
7 = 0.111 and (p * = 1.00, with Xini = 15.489 and </9 ini = 4.5. 


held model considered here. The sharp turn in the held trajectory occurring at 
N ~ 10 leads to the conversion of isocurvature modes into curvature modes, and 
essentially washes out any features in the scalar power spectrum. In turn, this 
results in the enhancement of the scalar power spectrum at k < 10 _2 Mpc _1 , 


leaving the horizon before the turn. It is important to note that the conversion 
of isocurvature modes in two-held models causes the curvature perturbations to 


be frame dependent {24~lj . We also see from Figure 3.11 that the tensor-to-scalar 
ratio over wavenumber k has the same form as that in a minimally-coupled single 


held case. 
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Figure 3.11: Power spectra (left) of scalar (V^, solid black) and tensor pertur¬ 
bations (Vt, blue dashed) and the associated tensor-to-scalar ratio r (right) 
against wavenumber k (Mpc -1 ) for the two-held model. The model parameters 
are m x = 5.6 x 10 -6 , m v = 4.48 x 10“ 5 , (3 = 0.009, 7 = 0.111 and (p* = 1.00, 
with Xi n i = 4.5 and <p- m \ = 15.489. 


3.3 Conclusion 

In this work, we have studied the effects of sharp transitions within the effective 
Planck mass during the inhationary epoch, specifically examining its effects on 
the primordial power spectra. We have focused on single-held models before ex¬ 
tending into a two-held model, which contains a secondary scalar held behaving 
as the inhaton held. In the single-held models, we find that such a transition 
gives rise to strong features in the scalar power spectrum, with these features 
occurring at the scales corresponding to those leaving the horizon during the 
transition, which are within our observational window. We see these features 
are rehected in the tensor-to-scalar ratio at the same scales. In contrast to 0 . 
we do not find oscillatory behaviour in the scalar power spectrum for these type 
of transitions. However this is to be expected as we do not study the same 
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theory. 

We have presented a single-field case in which both the Planck and BICEP2 
measurements agree, hence reducing tension between the two. Moreover, we 
have also studied the single-field case with the inclusion of a minimally-coupled 
field. In the two-field case, due to the conversion of isocurvature modes into 
curvature modes, we see dampened features in the scalar power spectra and 
in the tensor-to-scalar ratio. To further this study the generation of the CMB 
angular power spectrum and comparison to data will be required. 



Chapter 4 


Disformal couplings 


In the last chapter we studied scalar-tensor theories in the context of inflation; 
these theories can also be studied in the realm of dark energy. We will now 
make the transition to later times in cosmic history and consider the era of 
dark energy. Additionally in this chapter we will consider the background only 
and therefore, no perturbations will be discussed. 

As mentioned in Chapter 3, in scalar-tensor theories, there are two representa¬ 
tions of the action, namely the Jordan and Einstein frames. We now remind the 
reader of the frame definitions: the Jordan frame is a frame in which the scalar 
field and the Ricci scalar are coupled, whereas the Einstein frame is one where 
the Einstein-Hilbert action is as standard. Due to differences in representation, 
the question of which frame, Einstein or Jordan, is “most physical” arises; this 
question is addressed in the following references: [2441 245 . 1246 , . I247j . When 
presented initially with the Jordan frame, in order to obtain the Einstein frame 
the standard practice is to perform a conformal transformation of the metric; 
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conformal transformations have been studied extensively in [1581 11591124811249] . 
However, one can study a more general metric transformation between the two 
frames. 


It was first demonstrated in 1993 by Bekenstein [ 250 ] that the most general 
metric transformation between two frames which obeys both causality and the 
weak equivalence principle is 


g f w = C{4>)gnv + D{4 >)(. 


(4.1) 


This is known as the “disformal” transformation, which contains derivatives of 


the scalar field. The first term in Eq. (4.1) is the conformal coupling between 


the scalar and matter fields, and the disformal coupling is represented by the 
second term. Disformal couplings have attracted a significant amount of atten¬ 
tion since Bekenstein’s paper in 1993, [251112521125311254112551125611257] . These 
couplings have been studied as an alternative to dark matter, for example the 
Tensor-Vector-Scalar gravitational theory (TeVeS) [258]. It has led to the pro¬ 
posal of a new relativistic formulation of the modified Newtonian Dynamics 
(MOND) labeled bimetric MOND (BIMOND) [1259] . Other topics studied with 
disformal couplings are in bimetric theories of gravity with the inclusion of the 
Eddington-Born-Infeld (EBI) action [ 2601 261. 1262] . The inflationary scenario 
with the disformal relation has also been investigated [251] . The disformal re¬ 
lation has also appeared in variable speed of light theories [263 ;, [264]. 


In this chapter, we consider a scalar-tensor theory with disformally coupled 
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fields to study the effects of this type of coupling on the properties of the Cos¬ 
mic Microwave Background (CMB). In order to check consistency throughout 
this theory, both a fluid description and a kinetic theory was first established 
before discussing observables. Observables have been chosen to be expressed in 
terms of the metric in which an observer experiences, in this case, the baryonic 
metric. The baryonic metric is a metric in which baryonic matter is uncoupled 
to the scalar field. This choice was decided upon to avoid the frame-independent 
formulation for disformally coupled fields. The subject of frame-independence 
in conformally coupled theories is discussed in [265, [266 ]. 


The chapter is laid out as follows: in Section 4.1 the action and definition of a 
disformal coupling are presented, to be followed by the background equations 
of motion. Following this, the fluid description and kinetic theory are presented 


in Sections 4.2 and 4.3 A discussion on observable quantities is contained in 
Section [4.4[ leading into the presentation of results in Section 45 Finally, our 


findings are summarised in Section 4.6 


4.1 The model 


The action that we consider is 


S = I d A x ^g 


/If 2 1 


+ ^2 Si [g^jXi] > (4-2) 


where R is the Ricci scalar calculated with respect to the metric g^, Si denotes 
the action of the various matter fields Xi propagating on different geodesics via 
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the relation 

9%, = Ci((p)g > + Di^d^dvcj) . (4.3) 

It is important to note that this action is expressed in the Einstein frame. From 
now on, the reduced Planck mass are set to one, Mpi = 1. 


A cosmological setting is considered and the metric g [LV is of FRW form; for 


the FRW line element, we refer the reader to Eq. (1.8) in Section 1.3. By the 


variation of the action in Eq. (4.2) with respect to the scalar field and using 


Eq. (1.8) yields the Klein-Gordon equation for the held (j) 


4> + 3H<p + V' = , 


(4.4) 


and with respect to the metric g^ u , the Einstein field equations 




rjifU/ 


(4.5) 


with G^ v as the Einstein tensor and T^ representing the total energy-momentum 
tensor of the system. It is composed of two energy-momentum tensors: the 
scalar field and of the matter species 



rj-ifJ.U 


2 SS* 
V~9 s g^ ’ 

2 5Si 

V=gSg»» ’ 


(4.6) 

(4.7) 
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where S$ is the action of the scalar held. We shall assume that the energy- 
momentum tensor for each individual species is of perfect fluid form and hence 
T l ' u = diag(— p,p,p,p), where p and p are the energy density and pressure of 


the fluid respectively. Using this and the Einstein field equations in Eq. (4.5) 


in a FRW spacetime yields the conservation equation for the matter species 


Pi + 3 H(pi +pi) = -Qi<t >, 


(4.8) 


where the couplings, Qi are found by varying the action Eq. (4.2) are given by 


Qi = + |U, 4>„ rr - (§b,7T ] . 


(4.9) 


with T % representing the trace of the energy-momentum tensor of the species 
denoted by i. The following notation is used: dots represent derivatives with 
respect to cosmic time t, primes are those to the scalar field 0 and commas 
denote partial derivatives. The Hubble parameter is H — a/a. In the following 
we will consider two fluids labelled i — 1,2 and each is disformally coupled with 


a different magnitude given by D{. From Eqs. (4.8) and (4.9), the couplings Q\ 
and Q 2 are 


Ao 


q = ^ _ ( B B 2 D\Pi \ 

^ a,a 2 - d x d 2PxP2 V 1 a 2 ) ’ 

Q = A _ (b BiD2 PA 

A,A 2 - D 1 D 2 p 1 p 2 { 2 A, ) 


(4.10a) 

(4.10b) 
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where 


A.i 

Bi 


Ci + Di(pi — 0 2 ) , 

— ^ — 1 + 3—— —(f 1 
. 2 \ Pi) * 


DA3H 


R)< 


-)^ + v , + ( A ^ 2 


Pi 


(4.11a) 


(4.11b) 


In order to achieve the dynamics needed for slow-roll, a slowly varying potential 
is required. For this purpose, the exponential potential was chosen 


V(4>) = V 0 e- X * , 


(4.12) 


where Vq is a constant and A = 1. We now choose to study the effects of the 
disformal coupling without the presence of the conformal term; C(<f)) = 1, and 
D(y<p) is treated as a constant energy scale: = M” 4 where M is a constant. 


Before proceeding into the study of the CMB, two descriptions must be made 
consistent; these are the fluid description (used so far) and the kinetic theory. 
By doing this, we can recover the equations found through the variation of 


the action in Eq. (4.2) and have a more complete description of this disformal 
theory. 


4.2 Fluid description 

For simplicity, one fluid will be considered in this section; this can be easily gen¬ 
eralised to multiple fluids. We will revert back to the usage of multiple fluids 
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after this discussion. 


The fluid considered will have an arbitrary equation of state w and the disfor- 
mal coupling between the fluid and the scalar field is given by M, the constant 
energy scale. In the Jordan frame, the scalar field is decoupled from the matter 
fields and so the fluid considered is of perfect form 


V av = (p + p)u^u u + pg^ , 


(4.13) 


resulting in = diag(— p,p,p,p). The 4-velocity vectors of a comoving ob¬ 
server in the Jordan frame u M are defined as 


IF = 


dx^ 
df 



(4.14) 


where the relation between the proper time of the observer r and the time 
variable of the Einstein frame t is 


dr = 



(4.15) 


Likewise, the fluid considered in the Einstein frame is also of perfect fluid form 


= (p + p)F l u u + pg^ , 


(4.16) 
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with T IJ U = diiig(—p,p,p,p) and the Einstein frame 4-velocity vectors of a co¬ 
moving observer u M are 


^ = ^ = (1,0,0,0) . (4.17) 

Tildes will be employed to distinguish the Jordan frame from the Einstein frame. 
In addition to this, the contraction of metrics with various ranked tensors must 
be clarified. In the Einstein frame, all tensorial quantities are contracted with 
the Einstein frame metric 

A„ = g IJa A a , T^=g ua T^. (4.18) 

Similarly, all tensorial quantities in the Jordan frame are contracted with the 
Jordan frame metric 

A,= ~ gfia A a , t^=~g va t^. (4.19) 

The fluid in the Jordan frame satisfies the standard energy conservation equa¬ 
tion 

V /t T /i!y = 0 , (4.20) 

for which the covariant derivative is compatible with the metric g^, i.e. V^g^ = 
0. The energy-momentum tensor can be expressed in terms of the Jordan frame 
action: 

~ 2 JS 

rjifilP _ 

V-9 


(4.21) 



4.2. FLUID DESCRIPTION 


117 


At the level of the background, the scalar field is is only time dependent. 


A relationship between the Jordan and Einstein frame energy-momentum ten¬ 


sors can be obtained by using Eqs. (4.7) and (4.21) 


-g bguv 


= C\f l+ , 


(4.22a) 

(4.22b) 


and in the FRW spacetime, this relation is simplified to 


r r^ _ Wi _ y 

M 4 


(4.23) 


Note here that the dot represents the derivative with respect to t where goo = 
— 1, i.e. the time derivative in the Einstein frame. 


By simply taking the trace of the energy-momentum tensors for both frames 


and using Eq. (4.22b), the relationship between the two equations of state is 


P _ P, i 
P P 


M 4 


(4.24) 


The conservation equation in Eq. (4.20) can be rewritten using the relation 


between the two equations of state in Eq. (4.24), this yields 


p + 3 H I p + 


M 4 


M 4 - P ) M 4 - k > 2 


P 


(4.25) 
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Using the Klein-Gordon equation in Eq. (4.4) and Eq. (4.25), the conserva¬ 


tion equation found from using the variational principle stated in Eq. (4.8) is 
retrieved and hence, a consistent macroscopic description is achieved. 


4.3 Kinetic theory 

To make this disformal theory fully coherent, a microscopic description must 
be attained. In order to do this, a statistical approach is needed — see [ 2671 

12681125511270 ] . 

The energy-momentum tensor can be expressed in terms of an integral of the 
distribution function / over phase space. The distribution functions for the 
two frames, / and /, are given in terms of their position and momenta vectors 

[ 25511271 ] : 

dN = dx 1 dx 2 dx 3 dP 1 dP 2 dP 3 f (4.26) 

and 

dN = dx 1 dx 2 dx 3 dP 1 dP 2 dP 3 f , (4.27) 

where dlV and dlV are the number of particles in a differential phase space in 
the Jordan and Einstein frames respectively, dx 1 dx 2 da; 3 is the position phase 
space and dPidP 2 dP 3 is the momenta phase space for the Jordan frame (and 
similarly for the Einstein frame but with tildes). This formalism was also used 
for the conformal transformation case in 1272; . 


Likewise with the fluid description given in the previous section, each frame 
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has its own definition in kinetic theory (267, 268J : 



f d 3 P P^P U 

(4.28) 

rpiil> _ 

' V=9 P° 1 ' 


f d 3 P P^P V ~ 

(4.29) 

rpnv _ 

1 P° 1 ' 


Notice that the integrals are given in terms of the 4-momenta of that frame 
(tildes for Jordan and none for Einstein) and are as follows 




dx^ 


p/' 


dx 11 

dA~ 


(4.30) 


resulting in 


pM 


dx M dA dA 

_—— JDH _ 

dA dA dA ’ 


(4.31) 


where A and A are the affine parameters in the Einstein and Jordan frames, 
respectively. Within the Jordan frame, after decoupling the photons will obey 
the collisionless Boltzmann equation [2T8j 


d[ 

d t 


= 0 . 


(4.32) 


By using the relation between the two energy-momentum tensors Eq. (4.22a), 


their definitions: Eqs. (4.28) and (4.29) yields 



dA 

dA 


-9 


(4.33) 
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This relation can be applied to the distribution functions of the two frames 


Eqs. (4.26) and (4.27), and results in the following: 


d * = fxI dN - 


(4.34) 


From this relation, it is clear that the two quantities dlV and dlV are not frame- 
invariant and ultimately, in the Einstein frame the Boltzmann equation takes 
the form 


d 
d t 



0 . 


(4.35) 


The frame variance is due to the presence of the disformal coupling and is unlike 
to conformal transformations [265| . 


We now revert back to understanding how disformal couplings affects the prop¬ 
erties of the CMB; we know that photons will travel on geodesics in the Jordan 
frame and hence the following equations are true 

P^^P U = 0 , ~g^P v = 0 . (4.36) 


Note that the second of these equations can be rewritten as g flu P fi P l/ = 0 by 


using Eq. (4.31). Adding to this, we require the relation between the affine 
parameters defined in the Jordan and Einstein frame. In order to calculate 


this, we can use the fact that the energy-momentum tensors in Eqs. (4.28) and 








4.3. KINETIC THEORY 


121 


(4.29) can be rewritten for a system of point particles [274] 

pfi pv 


ptiu = 


P o 

pti pv 

P° 


J^(x-x') , 

(4.37) 

(S (3) (x-x') . 

(4.38) 


Similarly to the expressions of the distribution functions stated in Eqs. (4.26) 


and (4.27), the position vectors of the points x, are the same in both frames, 


and hence are the same Dirac delta functions. Using the energy-momentum 
tensor relationship between the two frames, their respective 4-momenta and the 


discrete point particle representation given in Eqs. (4.22b), (4.31), and (4.38), 


we find the relation between the Jordan and Einstein frame affine parameters 
to be 


rjnfiH _ 


1 _ 

M 4 




M 4 po 




rpfiv _ 


M 4 VdV ’ 


dA 

dA 


- = \1- 


M 4 


(4.39) 


To obtain various conservations through kinetic theory, the Liouville equation 
is needed 


cf = m , 


( 4 . 40 ) 
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where C is the Liouville operator and C[f] contains the collision terms. The 
Liouvillc operator acting upon a distribution function / is defined as 


d/ dx M df dP M df 
dA dA dxv dA dP^ 


(4.41) 


The geodesic equation in Eq. (4.36) can be rewritten as 


d P^ 

7T 


+ f 


0, 


(4.42) 


where the connection is given by the Christoffel symbols. In the equation given 
above, Ffp are the Christoffel symbols in the Jordan frame. For completeness, 
when C(f>) = 1 and D(f>) = M -4 the nonzero Jordan frame Christoffel symbols 
in terms of the Einstein frame are 


f-0 _ # 

00 M 4 — 0 2 ’ 

(4.43a) 

f°, = 1/1 . a 2 H5 

13 M 4 - 0 2 1 

, (4.43b) 

f} 0 = H5) . 

(4.43c) 

Using Eqs. (4.31), (4.41) and (4.42) — expressed 

in the Jordan frame, we can 


write down the Einstein frame Liouville operator: 


Cf 


d/ 

dA 


ai_ I poflfL) 

Or" ] dA 2 \dA/ 


, -nO 

^ 1 a/3 



(4.44) 
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Now referring back to the Einstein frame Boltzmann equation in Eq. (4.35), we 
can use this to calculate (d//dA) 


0 


d/ 

d t 
df_ 
dA 


/dA 

dt VdA ; 


d 

dt 


dA 

dA 


'dA 
' dA 


dt 


'dA 
' dA 


dt 


/*£(£)■ 


-/• 

(4.45) 

~P°f ■ 

(4.46) 


Applying the relation between the two affine parameters in Eq. (4.39) and the 


Christoffcl symbols stated earlier to the Liouville equation in Eq. (4.40), yields 
the Einstein frame Boltzmann equation 


P° 


d _l 

dt 


— HPP° 


d 

dP 




(4.47) 


where we have also used the geodesic equation g flu P^P l/ 


0 Q which leads to 


d P° _ M 4 a 2 P 

d P ~ M 4 - 0 2 P° ’ 


(4.48) 


with P 2 = 8ijP l P j . We will require definitions of the energy density, pressure 
and number density derived from the statistical representation of the energy- 


remind the reader that we rewrite the geodesic equation in Eq. (4.361 by using 
Eq. Oil 
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momentum tensor; they are given as follows 


P = 

P = 

n = 



- P°f 

d 3 P a 2 5 ij P i P j 

J ’ 

d 3 P 

y/=9* ’ 


(4.49a) 

(4.49b) 

(4.49c) 


where n is the particle number density. By integrating Eq. (4.47) in 


space, we obtain again the same equation as Eq. (4.25) and hence, 
theory is consistent with both the fluid and action approaches. 


momentum 
this kinetic 


When Eq. (4.47) is compared to the the standard collisionless Boltzmann equa¬ 
tion, we see an extra term present on the right-hand side. The presence of this 
additional term is a result of photons not travelling on geodesics in the Einstein 
frame. We can interpret this modification as an “effective” collision term as it 
contains no derivatives of the distribution function /. 


In addition to its ability to produce the conservation equation, the Liouvillc 


equation in Eq. (4.40) can also be used to construct an equation describing the 


evolution of the particle number density 


[ d 3 P tf f d 3 P C[f] 

J V=gP° J V~9 p° 


(4.50) 


with attention placed upon the factor of 1/P° on both sides of the equation. 


Applying this factor to the Einstein frame Boltzmann equation in Eq. (4.47) 
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yields 


— — H pR- — 

dt dp 


M 4 -(p 2 


/, 


(4.51) 


and by integrating Eq. (4.51) over momentum space whilst utilising the definitions 


stated earlier in Eqs. (4.49) yields the equation for the particle number density 


in the Einstein frame, which is as follows 


h + 3 Hn = > 


(4.52) 


with F = ( d\/d \) 2 . This equation implies that the particle number in the 
Einstein frame is not conserved and hence the CMB is not a blackbody (and 
does not obey the adiabaticity condition in Eq. (19) in [275] — this relation is 
also used in m\)- As a further check, the equation of state for the Einstein 
frame obtained from kinetic theory is 


9^T 


av 


—p + 3 p 
-p + 3 p 
3 p 


9nc 


9iio 


^3p p)ctjph- 

V^g p0 


'ifi 'ria 


/, 


M 4 

d 3 P 


d.3 p p a p v 

V^g p0 


/, 


y/^g 


P°f , 


4> 



(4.53) 


where g fiU P f ^ l P u = 0. We find it is the same as that from the fluid description 


Eq. (4.24). 
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From this, we can state that the equation of state is not frame-independent in 
the presence of disformal couplings. This is in exception to when the pressure in 
the Jordan frame is zero, which in turn makes the Einstein frame pressure zero. 
In addition, the distribution function in the context of disformal couplings, is 
also revealed not to be a frame-independent quantity. 


4.4 Observables 

We shall return to the original problem with two species: radiation and matter, 
with matter containing both baryonic and dark, that are separately disformally 
coupled to the scalar field. 


In this scenario, it is important to note that there are three distinctive frames: 


The Einstein frame, which takes the action of the form stated in Eq. (4.2). 


• The radiation frame. This is a frame in which the CMB radiation is un¬ 
coupled from the scalar field and due to this, matter is in general coupled 
to the field. 


• The Jordan frame. In this frame, radiation is in general coupled to the 
scalar, with matter (encompassing both baryonic and dark) uncoupled. 


For the final two frames described, the gravity-scalar part of the action will have 
a nonstandard form. Note that if the strength of the disformal couplings in both 
matter and radiation are the same, both the radiation and Jordan frames co¬ 
incide. From the previous section, we know that the distribution function and 
the equation of state is different for each frame, and this affects observables. 
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We will perform all calculations in the Einstein frame for simplicity. Once this 
is done, observables will be expressed in terms of the Jordan frame in order to 
compare with data. The Jordan frame is the most natural choice as this is the 
frame where matter (and hence experimental equipment) is uncoupled. 


To search for disformal couplings we employ two tests. Two such tests cho¬ 
sen are the redshift-temperature and the /./-distortion; both the temperature of 
the CMB radiation and spectral distortions have been previously used to place 
constrains on extensions to the standard theory [ 276 . [277]. 


The CMB temperature evolution has been studied in the redshift range 0 < 
z < 3 by the following authors [UEJE, 278]. In order to study the temperature- 
redshift relation, an expression of the redshift is required. The definition of the 
redshift is given in terms of the fluid’s velocity vectors and 4-momenta; it is 
defined as m 


l + z = 


(■ U nP 1 *) 


(m) 

em 




(4.54) 


where z is the redshift, u^ are the 4-velocity vectors and are the 4-momenta. 
The subscripts “em” and “obs” stand for emission and observation respec¬ 
tively, with superscript “m” representing the matter frame. In this expression, 
all quantities are those which belong to the Jordan frame; the 4-velocity is 
u /t = (—c 0 bs, 0, 0, 0) and P p is the measured 4-momenta of the photon. The 


measured speed of light c G b s is given in Eq. (4.64). 


Using the requirements of radiation in the Jordan frame (Eq. (4.36)), we know 
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that photons will follow disformal geodesics. Equation (4.36) in combination 


with Eq. (4.54) reveals that the redshift relation in terms of the baryonic metric 


is the same as that found for General Relativity: 


l + z = — , (4.55) 

a 

where a 0 is the present day scale factor. An effective temperature is assigned 
to the energy density of the CMB radiation in the Jordan frame p y m \ through 
the relation p oc T 4 (279]. 


The second test involves the study of ^-distortions. Previously, we stated that 


the distribution function / alone is not conserved but instead obeys Eq. (4.35). 
We make the assumption that the CMB when first produced is to a very high 
accuracy a blackbody. The distribution function in the radiation-dominated 
epoch is of Planckian form 


/K r > = ■ < 456 ) 

and in addition to this, 0 = 0. We further ignore the, if any, spectral distortions 
that were generated during the early Universe. Over time the contribution made 
by the derivative of 0 will become more significant and will affect the form of 
the distribution function: 

/obs=/im Gi) ’ (457) 

where f Q b s and / ini are the observed and initial-time distribution functions and 
subscript “0” indicates at observation i.e. at present time. The presence of 
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the factor dA/dA can be quantified as an effective chemical potential g in a 
distribution function. Hence, the observed distribution function can be written 
in the form 

/K T ) = e ,/rt _ ! ' ( 4 -5S) 

The chemical potential should be viewed as a way to quantify the deviation 
from the Planckian spectrum caused by the disformal coupling. 


Given that observationally ^-distortions are very small i.e. g -C 1 — obser¬ 
vations include those made by COBE FIR AS [7[ 07] and PIXIE [ 280] . see Dent 
et al. for a discussion on /i-distortions BED - we can use this in combination 


with Eq. (4.35) to find that 


g ~ (c obs - 1)(1 - e u/T ) , 


(4.59) 


where c D b s is the measured (dimensionless) speed of light. By using the following 
line elements for the observer and photons 


dsL = 9™vdx^dx" = -d t 2 + dx 0 2 bs , (4.60) 

ds 2 = g^dx^dx" = 0 , (4.61) 


and the metrics as stated in Eq. (4.3) with C — 1 and Di = M i where 


i indicates the matter species i.e. matter and radiation, we can rewrite the 
photon metric in terms of the matter metric 
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From the line elements stated in Eqs. (4.60) and (4.61), this leads to 


ds 2 = 


1 flV 


M 4 M 4 

7 m 


d ^obs 

d t 2 

V-l I y,-, 


-"obs 


= 1 ~ 


= 1 - 


+ da? 0 bs + [ 


1 

1 N 

\# 2 

M 4 


/ dr 2 ’ 

1 

1 ' 


M 4 


1 Vd r m 


i 


dx^dx v 


d 0 2 , 


(4.63) 


where we have used the fact that when an observer is at rest ds 2 bs = —dt 2 = 
—d t^, and r m is the proper time of the observer. Alternatively, the observed 
speed of light can be expressed in terms of Einstein time variables 


1 - 0 2 /M 7 4 

i - y/M,t' 


(4.64) 


There are limits which have been experimentally found by COBE and FIRAS 
[7]: |/x| < 9 x 10~ 5 . This will be used to constrain the M m x M 7 parameter 
space. 


4.5 Results 


The Klein-Gordon equation and the conservation equations for the various 


species in Eqs. (4.4) and (4.8) are integrated such that the present day density 


parameter values are fl m) o = 0.3 (cold dark matter and baryons), = 0.7 


(dark energy) [92j and current temperature of the CMB is T 0 = 2.725K [4j. 
The initial field value is fixed at = 1.5 and Vq is allowed to vary in order to 
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obtain the our present day Universe. The two methods stated in the previous 
section: temperature-redshift and //-distortions, are then calculated. 


The expected temperature evolution as a function of redshift is numerically 


calculated for a range of M 7 ; our results are shown in Figure 4.1 Figure 4.1 
displays the values of M 1 which are in agreement and those not with current 


measurements from 0 0 El EZE!- Also included is the prediction of General 
Relativity, for M —y oo. A large M 7 results in a weak coupling between the 
scalar field and radiation. We find that this weak coupling distorts the linear 
temperature-redshift relation found in General Relativity. However, the cou¬ 
pling has a damping effect on the field which increases as the coupling increases 
(equivalent to M 1 decreasing). When the coupling is very large, the dampening 
is so severe that it results in <fi —>■ 0. Due to disformal couplings being composed 
of field derivatives, it is clear that when this limit is reached the coupling will 
vanish (shown by the tripledot-dashed line in Figure [4T|. 


A chi-squared minimisation method with 68% (la) confidence level is used 
to obtain an excluded region in M 7 x M m parameter space. By comparing the 


expression for the //-distortion in Eq. (4.59) against measured constraints in [7], 
we find them to be similar when under the assumption that the frequencies v 
are of order T 0 . Bounds on the M m x M 1 parameter space acquired from both 


methods are presented in Figure 4.2 


We find that the constraints found from //-distortions are much stronger than 
those from the temperature-redshift relation. Adding to this, the overall shape 
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Figure 4.1: Plot of the ratio T(z)/T 0 where T 0 = 2.725K, against redshift, z, 
for the exponential potential. The black circle marked at z = 0 is measured 
by COBE [4j. Different sets of measurements and their corresponding errors: 
Those marked with crosses are given by [5] and triangles are by |6]. There 
are five lines marked on this plot and each has an associated M 7 : Solid line 
is M 7 —» oo, dashed line M 7 = 2.203 x 10” 5 eV, dashed line with one dot 
M 7 = 3 x 10~ 5 eV, dashed line with two dots M 1 = 1.5 x 10~ 3 eV and dashed 
line with three dots M 1 = 2.2188 x 10~ 3 eV. For this plot, we fix M m = 0.05 eV. 
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log (M m / meV) 


Figure 4.2: Bounds on the M m x M 7 parameter space. The solid diagonal 
line represents M m = M 1 . The darkest shaded region is excluded by the CMB 
temperature evolution alone, we show the exclusion region above 68% C.L. . 
The light-gray shaded regions are excluded by the measured constraints on the 
/i-distortions [Tj. The line-shaded (hatched) area is excluded from our search 
of the parameter space as the numerics are unreliable in that region. Note that 
the regions excluded by the /x-distortion never touch the line M 7 = M m , since 
in this limit the /i-distortion vanishes. 
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from both methods are alike, which draws the conclusion that there must be a 
difference in the two couplings for any disformal effect to be observable. 


4.6 Conclusion 

From this study of disformal transformations, there are two points worth em¬ 
phasising: firstly, we find that for a given fluid, the equation of state under dis¬ 
formal transformations is in general, a frame dependent quantity and secondly, 
this coupling has an effect upon the distribution function, and consequently 
making the distribution frame dependent. By considering CMB radiation con- 
pled disformally to a scalar field, we found a modification to the distribution 
function of the CMB photons. This modification disappears when both baryon 
and radiation coupling strengths to the scalar field, are equal. 

Throughout this work, all observable quantities were expressed in the Jordan 
frame; in this frame, the baryonic matter and scalar field are uncoupled. This 
work can be further developed. For instance, the couplings strengths for both 
M m and M 1 should be in general, allowed to be dependent on the scalar field. 
At the start of this work the conformal factor was set as C(4>) = 1, with the pur¬ 
pose to study a purely disformal scenario. The most natural follow-up to this 
study is to examine the scalar field dynamics and ultimately, the constraints 
posed by the temperature-redshift relation and //-distortion, with the inclusion 
of the conformal coupling. In addition, conformal couplings require screening 
in higher density regions of the Universe clue to a fifth force arising in such 
models, so naturally disformal couplings with a screening mechanism should 
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also be studied [256] . In these models, the quantity cj) is expected to be much 
smaller and this will in turn, diminish contributions from the disformal cou¬ 
pling. A study into the evolution of the first order perturbations is essential as 
disformal couplings will affect the CMB temperature anisotropies and structure 
formation. 
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Chapter 5 


Conclusion 


There are two periods of accelerated expansion within the lifetime of our Uni¬ 
verse: the first is inflation which occurs at early times, and second is a result 
of dark energy at later times. Currently standing, there is an abundance of 
concrete models that are able to describe the dynamics during both of these 
periods. Many theoretical cosmologists have developed theories including those 
from string theory. One type of theory considered to be a more natural ex¬ 
tension of General Relativity are scalar-tensor theories. In this Thesis we have 
three cosmological studies with the first two in the context of inflation, and the 
final one in a dark energy scenario. 

In Chapter 2, three models of inflation were considered with all containing 
at least one nonminimally coupled scalar field. The models taken into account 
included double inflation, followed by one with a kinetic coupling and lastly, 
a DBI field; these were chosen with the purpose to study the effects of their 
kinetic terms upon the power spectra evolution for the curvature and isocurva- 
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ture perturbations. We found that the inclusion of these noncanonical kinetic 
terms causes the isocurvature power spectrum to vary. To be specific, if infla¬ 
tion is driven by the scalar held containing such noncanonical kinetic term, this 
would cause deviations of the isocurvature power spectrum amplitude from the 
case of double inflation. From our study, we found that a DBI held driving the 
inhationary period gives rise to the largest isocurvature perturbations, which 
will in turn affect predictions of the B-mode polarization in the CMB, which are 
currently being under investigation by experiments such as BICEP2 and Planck. 

In the following chapter, we studied the scalar and tensor power spectra gen¬ 
erated during inflation with the effects of a sharp transition in the effective 
Planck mass. We performed all calculations within the Jordan frame, i.e. phys¬ 
ical frame. In the single-held model we considered, the effects of the sudden 
transition can clearly be seen on the scalar power spectrum, however, not so 
evident on the tensor spectrum. We have also shown that with the inclusion of 
this transition occurring at scales close to those leaving our horizon, and hence 
within our observable window of scales, it is possible to reduce tension between 
the Planck and BICEP2 observations of the tensor-to-scalar ratio. To further 
this work, the addition of a minimally-coupled scalar held was made. However, 
the presence of the second held caused the suppression of all features resulting 
from the sharp transition, within the primordial power spectra. 

In Chapter 4, we move forward to later times in order to study dark energy 
with the addition of a nonconformal coupling, namely the disformal coupling. 
We considered the model in which all matter species, including radiation and 
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matter, are disformally coupled to the scalar field. To ensure the theory is con¬ 
sistent throughout before pursuing observable quantities, the fluid description 
and kinetic theory were established. From this, we found that the equation of 
state for the species in question, is a frame dependent quantity. In addition 
to this, the distribution function of the CMB changes between frames under 
disformal transformations (but is invariant under conformal transformations). 
This modification ends when both the radiation and matter species are coupled 
with equal magnitude to the scalar field. In order to relate this coupling to 
observations, we choose to perform calculations in the physical frame. For the 
case where only disformal couplings are present, we have also placed constraints 
on the strengths of both the radiation and matter couplings through both the 
temperature-redshift relation and //-distortions. 

To extend the study in Chapter 2, the full primordial power spectra for the 
curvature and isocurvature modes for all models considered will need to be 
produced, and this to be followed by research into the pre- and reheating mech¬ 
anism. Post-inflation dynamics will also be required for the theory studied in 
Chapter 3. One logical extension to the project on disformal couplings, in par¬ 
ticular to radiation, is to consider the perturbations generated. All the models 
considered in this Thesis will require the generation of the CMB angular power 
spectrum, and then for this to be compared to experimental data. 

Within the last three years, there have been many experiments releasing their 
data including those from Planck and BICEP2, all of which will place further 
constraints on theoretical models, and possibly rule some out. The results of the 
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BICEP2 experiment state that B-mode polarization has been detected. How¬ 
ever, there are questions about their analysis and as such, has led to suspicions 
about their findings on gravitational waves. At this current time, Planck has 
not released their results on polarization, but are in slight disagreement with 
those from BICEP2. With the first two inflationary studies presented in this 
Thesis, we see that it is possible to find models producing isocurvature modes 
relating to these B-mode polarizations, and it is possible to resolve disagree¬ 
ments between the two major experiments. For the case of dark energy, we still 
have no clue what it is, and with the proposal of unusual couplings between 
matter species, it will lead to more research within this area. With all this said, 
we will always require experiments with ever increasing precision to probe the 
CMB, and hopefully reveal that our Universe is not governed by the dynamics 
of minimally-coupled scalar fields. 



Appendix A 


Numerical method 


This appendix describes the numerical methods that were required throughout 
this Thesis. We have used the method outlined by [ 216] . but there are many 
other approaches available in the literature, for example see references [2821.199] . 


A.l Initial conditions 

The initial conditions for the background fields are chosen such that they pro¬ 
duce an inflationary period which lasts for a minimum of 60 e-foldings; the end 
of inflation is recorded when the slow-roll parameter e reaches 1. In addition, 
the derivatives of the background fields are given their slow-roll values. The 
Bunch-Davies vacuum fluctuation of the scalar field in terms of conformal time 
is given by {68] 

S>Pm = —Ue"'*" , (A.l) 

a\J2k 
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and is the oscillatory initial condition used in the numerical runs, the derivative 


of Eq. (A.l) with respect to conformal time is 


¥'ini = 


iV2k 


(PL + ik)e 


—ikrj 


(A. 2 ) 


where PL = a'/a is the conformal time Hubble parameter. For the case of a DBI 
held, this initial conditions changes slightly to include the Lorentz-like factor 7 

m 


<5Xini = — k=e ’i" , (A.3a) 

a7V2/c 

^ w= “^( K+ 7 + b) e "‘‘ v (A ' 3b) 


A.2 Set up 


Most, if not all differential equations mentioned in this Thesis e.g. background 
Klein-Gordon and held perturbation equations, are of second-order. I 11 order 
to solve these equations numerically, each second-order differential equation 
will need to be rewritten to form two hrst-order differential equations. The 
commonly used integration variable is the e-fold number as stated in Chapter 1 
where N = In a. For example, to solve the background Klein-Gordon equation 
for a minimally coupled scalar held with respect to the number of e-foldings, 
the variables are 
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and the set of differential equations which need to be solved are 


y ,n = 


/ 

V 


4>,N 


\ 



YM. 
h 2 y 


(A.5) 


where indicates derivative with respect to e-fold number. The differential 
equations in this Thesis were solved using Runge-Kutta 5th order and 8th or¬ 
der methods. It must be noted that for the case of the DBI held, it is best to 
perform the integral in logarithmic time, x = lnt. 


The code begins with the condition 


k 
a.H 


> 1 , 


(A.6) 


which ensures that the held perturbations are deep within the horizon. Note 
that the wavenumber k is the same as that arising in the perturbation equations. 


We treat the perturbations as stochastic variables deep within the horizon. In 
addition to this, the modes must be uncorrelated deep within the horizon and 
to simulate this we must perform two numerical runs [ 216] : the hrst run is set 
where the gauge-invariant perturbation of the held </> is set in the Bunch-Davies 
vacuum and the gauge-invariant perturbations of the held x is set to zero, and 
for the second run the conditions are swapped. 


Each of these runs are used to evaluate the curvature perturbation, either ( 
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or IZ, and then combined to form the curvature power spectrum 

u 3 

n = ^(ICi | 2 + K 2 | 2 ), (A.7) 

where the subscripts 1,2 denote the run number. This method is the same for 
the entropy power spectrum. 
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